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An introduction to the transport properties of graphene combining experimental results 
and theoretical analysis is presented. In the theoretical description simple intuitive 
models are used to illustrate important points on the transport properties of graphene. 
The concept of chirality, stemming from the massless Dirac nature of the low energy 
physics of the material, is shown to be instrumental in understanding its transport 
properties: the conductivity minimum, the electronic mobility, the effect of strain, the 
weak (anti-)localization, and the optical conductivity. 
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I. INTRODUCTION 

For a long time, theorists and experimentalists alike 
have considered the existence of a true two-dimensional 
(2D) material, having the thickness of a single atom - 



a one atom thick membrane - to be impossible. The 
reasoning behind this statement relies on the fact that 
both finite temperature and quantum fluctuations con- 
spire to destroy the otherwise perfect 2D structure of the 
hypothetic material. These fluctuations, originated from 
atomic vibrations perpendicular to the plane of the ma- 
terial, would preclude the existence of a true flat phase 
and concomitantly the existence of such a system. 




Figure 1 (Color online) An optical image of a graphene flake, 
obtained form the exfoliation of graphite, with an area of < 1 
mm^, on top of a silicon oxide wafer (courtesy of P. Blake). 

Nevertheless, in 2004, a group led by A. K. Geim, from 
the University of Manchester, U.K., isolated such a 2D 
material fNovoselov et al. 2004 2005b') . Under the name 



of graphene, this new material is an allotropic form of car- 
bon, with the atoms arranged in a 2D honeycomb lattice. 
The reason for the success lies on the isolation method. 
The developed method permitted one to isolate the 2D 
material on top of a 300 nm thick wafer of silicon oxide. 
The weak van der Waals interaction induces adhesion 
between graphene and the wafer, and once on top of the 
wafer, it is possible to move about the 2D material, trans- 
ferring it from one substrate to another, or even having it 



suspended over a trench, supported from one side (Booth 



et al. 2008 1. In the production method, graphite plays a 



key role, since this 3D material is itself made of stacked 



graphene planes (binded by van der Waals forces). The 
ingenuity of the method was then to find a way of peel- 



ing a single layer of graphene out of graphite ( Novoselov 
et al.\ |2004[ |2005b| ). Up to this date, the exfoliation of 
graphite can produce graphene crystallites as large as ~1 
mm^ (see Fig. [l|. The study of graphene became, since 
2004, an active field of research in condensed matter, 
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which holds many promises (Castro Neto 

Neto et aL| [20061 IFuhrer e^ aY| |2010| |GeimH2009't|Geim 

and Kim. 2008 Geim and MacDonald. 2007 Geim and 



Novoselov, |2007i |Katsnelson, ,2007, |Peres, ,2009, jService 
2009F^ 



Being the first truly 2D material, it is natural to ask 
how its properties differ from those of more conventional 
systems, such as the 2D electron gas in the inversion 
layer of an ordinary semiconductor. The current efforts in 
graphene research have focused on the interplay among 
elastic, thermal, chemical, and electronic properties of 
the material, with a special emphasis on charge and heat 
transport, and on optical properties. The need for a deep 
understanding of the transport properties of graphene 
is obvious, since the material is a potential candidate 
for incorporating the future generation of nanoelectronic 
and nanophotonic devices ([Blake efal. 2008,|Liao et al 



20T0[ [Lin e^ all [20T0l [Mueller et aZ.[[2010l[Schwierz[[2010 



Xia et al. 2009). Also in biophysics, graphene is finding 



new applications (Lu et al. 2010 Schneider et al. 2010 



Wu et al. 2010). Additionally, and of no less importance. 



graphene provides a realm for the emergence of new and 
exciting physics. 

In the field of electronic applications, faster electron- 
ics requires smaller devices, in particular because at the 
nanoscale it may be possible for the electrons to travel 
across some of the components of a device almost unim- 
peded. In a normal conductor, one of the sources of elec- 
trical resistance is scattering of electrons by impurities 
and defects (and at room temperature, also by phonons). 
A measure of the effect of impurities on the electronic 
transport is the mean free path i (the average distance 
traveled by an electron between two consecutive colli- 
sions), which in a material with high degree of purity 
and with small dimensions can be larger than the typical 
length of the system L^ leading, in these circumstances, 
to what is called ballistic transport (in this regime the 
current becomes spatially non-uniform) . It just happens 



that in graphene £ can be as large as 1 /j,m (Bolotin et al. 
2008a Novoselov et al. 2004| ), putting graphene into to 



the ballistic regime, since the typical size of graphene 
based field effect transistors is L^ ^0.25-0.5 /im (Du 



et al.\ 20081. 



The first ground breaking publications of the Manch- 



ester's group (Novoselov et al. 2004 2005b I not only 



made the method of isolating graphene immediately pub- 
lic, but also established the major relevant problems in 
graphene transport: the ambipolar field effect (see Fig. 
^ , the independence of the electronic mobility upon the 



gate voltage, the large electronic mean free path, the con- 
ductivity minimum and the absence of Anderson localiza- 



tion ( Bardarson et al. 2007 1 , the magneto-resistance and 



the chiral quantum Hall effect (Novoselov et al. 2005a 



Zhang et al. 2005). These topics still orient much of the 



research in graphene physics at present. 

Since the publication of a comprehensive review on the 



theoretical properties of graphene (Castro Neto et al. 



2009), there has been additional relevant contributions 



to experimental and theoretical studies of its transport 
properties. In this Colloquium, we present an update 
on the experimental and theoretical developments in this 
fast growing subfield of graphene research, at a level ap- 
propriate to graduate students entering the field. 



II. THE ENERGY SPECTRUM OF GRAPHENE AND THE 
EMERGENCE OF DIRAC ELECTRONS 

As stated, graphene is a 2D material made solely of 
carbon atoms, arranged in a hexagonal lattice such as 
that shown in Fig. ^ There are five vectors represented 
in Fig. |2j the three next-nearest neighbors vectors di 
{i = 1,2,3), and the primitive cell vectors ai and 02. 
We further note that the hexagonal lattice is made of two 
inter-penetrating triangular Bravais lattices. Therefore, 
the effective model describing the low-energy physics of 
graphene has to keep track of the two atoms per unit 
cell, characteristic of the honeycomb lattice. Electrons 
in graphene can be described by a tight-binding Hamil- 
tonian reading (spin index omitted) 



H=-tJ2\ARn){Rn+S„B\+E.C. 



(1) 



n^Si 



where |A, i?„) represents the Wannier state at the unit 
cell Rn- and the equivalent definition holds for jS, fl„ + 
5i)\ t is the hopping energy. This Hamiltonian describes 
the motion of electrons in the tt— orbitals of the material, 
made from the hybridization of the atomic 2pz orbitals, 
and includes both low-energy and high-energy electron 
states. The calculation of the electronic energy spec- 
trum of graphene proceeds by introducing, in Eq. (IT]), 
the Fourier representation of the Wannier states in terms 
of the Bloch states of momentum fc; the spectrum then 



reads ( ,Castro Neto et al.\ [2009[ [Wallace[ [1947[ ) 



Eik) =±t\\ + e 



ikai 



^ih-CL2 \ 



(2) 



It is immediately obvious that the band structure of the 
TT— electrons is composed of two bands, one at negative 
energies (a hole band) and the other at positive ones (a 
particle band). In the Brillouin zone there are two spe- 
cial, non-equivalent (i.e. not connected by a reciprocal 
lattice vector), wave numbers, termed K and K' , and 
shown in Fig. |2] The transport properties of graphene 
are mostly determined by the nature of the spectrum 





Figure 2 (Color online) Real space lattice and Brillouin zone 
of graphene. Left: The hexagonal lattice of graphene, with 
the nearest neighbor di and the primitive a^ vectors depicted. 
The area of the primitive cell is Ac = 3\/3cto/2 — 5.1 A'^, and 
ao — 1.4 A. Right: The Brillouin zone of graphene, with the 
Dirac points K and K' indicated. Close to these points, the 
dispersion of graphene is conical and the density of states is 
proportional to the absolute value of the energy. 



around these two points. Close to K and K' the disper- 
sion, Eq. (pi, is conical, and given by E{k) = ^.vpfik, 
with vp — 3taQ/2h, where k is the momentum mea- 
sured relatively to either K or K', depending on the 
position of the cone in the Brillouin zone. Using the 
widely accepted value of t ~ — 2.7 eV for the hopping 
(in reality the values of t vary in the literature, spanning 
the interval from -2.7 to -3.1 eV) we obtain vp < 10^ 
m/s. The experimental studies are consistent in obtain 
ing Vp ~ 1.1 X 10^ m/s (Jiang et al. 



et al. 2005a Zhang et al. 2008 20051. A direct mea 
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surement of the Dirac spectrum in graphene has recently 
been obtained using angle-resolved photo-emission spec- 
troscopy (Sprinkle et al. 20091. Since each carbon atom 
(electronic configuration Is^ 2s^ 2p^) hybridizes with its 
three nearest neighbors according to the hybrid orbitals 
sp^, there is one electron left in the pz orbital. Therefore, 
the system is half filled, with the important consequence 
that the low-energy physics is controlled by the spectrum 
close to the K and K' points. Many of the new and ex- 
citing properties of graphene stem from this fact. The 
vicinities of these two points are also referred to as the 
two valleys of the electronic spectrum of graphene. 

The spectrum E(k) — ±vphk is formally equivalent to 
that obtained from solving the 2D massless Dirac equa- 
tion. Indeed, it is easy to show (Castro Neto et al. 2009 



Hamiltonian for the electrons in graphene has the form 



H 



K = vpcr p, 



(3) 



whereas close to K', the Hamiltonian is obtained from 
Eq. (pl by making the transformation Hk' = —Hk- 
The operator cr is written in terms of the Pauli matrices 
as cr = (ux, (Jy), and p is the momentum operator. Com- 
puting the eigenvalues of the Hamiltonian (Is]) , the conical 
spectrum indicated above is immediately obtained. We 
stress that cr does not represent real electronic spin; it 
is instead a formal way of taking into account the two 



carbon atoms per unit cell in graphene, as we have an- 
ticipated above. For this reason, cr is termed pseudo- 
spin. The density of states associated with the coni- 
cal dispersion of electrons in graphene is computed by 
determining the number of states per unit cell in the 
Brillouin zone N{E) up to the momentum k. Taking 
into account contributions from states near K and K' 
points, we obtain N{E) — fc^Ac/(27r), from which the 
density of states p{E) per spin and per unit cell is given 
by p{E) = dN{E)/dE = 2\E\/{-K^t^) , and the prim- 
itive cell area, A^ is defined in the caption of Fig. |2] 
The linear dependence of the density of states on energy 
is one of the fingerprints of massless Dirac electrons. For 
neutral graphene, the Fermi energy is zero. Therefore, 
the density of states vanishes in this case. 

The electronic linear spectrum and the chiral nature 
of the electron's wave function (see below) make elec- 
tronic behavior in graphene quite unique, and are re- 
sponsible for the remarkable properties of this material. 
Since cr ■p\'ip) = ±p|'(/'), then the operator h — cr -pjp has 
only two eigenvalues ±1. The operator h is known as the 
helicity operator, and has the following physical inter- 
pretation: in an energy eigenstate, the pseudo-spin cr is 
either parallel or anti-parallel to the momentum p. In the 
K. valley, electrons have positive helicity and holes have 
negative helicity; in K.' the opposite happens. The he- 
licity (or chirality) of electrons in graphene is responsible 
for the Klein tunn eling effect (|Beenakker[ |2008[ |Cheianov| 
[md Fal'kol |2006| [Katsnelson et aZ.| |2006|), observed re 



cently in graphene heterojunctions (Stander et al. 2009 



Young and Kim 2009[). We then see (and at odds to high- 



energy neutrino physics) that massless Dirac electrons in 
graphene come with both right and left chirality: parity 
is a symmetry of graphene. In Fig. [3] we show, in sim- 
ple terms, the origin of the Klein tunneling effect: the 
probability of electronic transmission through a poten- 
tial barrier is equal to 1, for head-on collisions; it is said 
that backscattering is suppressed. 

We should note that chirality is not, however, an exact 
symmetry of the problem. This occurs because the spec- 
trum of graphene is not exactly linear at all energies. 
The deviation from the perfect massless Dirac behav- 



ior is known as trigonal warping (McCann et al. 2006 



Semenoff 19841 that close to the K point the effective £■ > 1 eV 



Narozhny 20071, and starts playing a role for energies 



We remark, however, that trigonal warping 
might be important for observation of weak localization 
at energies much lower than 1 eV (see Sec. V.Bl. 

The solution of the eigenproblem iJjfl^) = E\i\}) is 
easily obtained by recognizing its formal equivalence to 



that of a real spin in a magnetic field ( Castro Neto et al. 



2009| , with the wave function reading 

g-j0(fe)/2 



|V'±) = 



1 



r 



ih-r 



M±(fc)e"' 



(4) 



72 V ie^'^'^'Z" 

and 0(k) = aictan^ky / k^) . Since the eigenproblem we 
have just solved is formally identical to a spin one-half 




Figure 3 (Color online) At the K valley, electrons have pos- 
itive helicity, h = 1, whereas at the K' one, the helicity is 
negative (Q = K' — K represents the transferred momentum 
when a scattering event between the valleys takes place). In 
a head-on collision of the electron on a potential barrier, the 
backscattered electron has to change its momentum from q 
to —q. For such a head-on collision (taken here along the 
a;— direction), ^ is a constant of motion, with eigenvalue -1-1, 
but backscattering would imply a modification of this eigen- 
value to —1. This, however, cannot be, because h is a. con- 
served quantity, then the transmission probability through 
the barrier, for such type of collision, has to be one. Thus, 
backscattering is suppressed for intra- valley scattering events. 
On the other hand, electrons in the K' and K valleys have 
opposite chirality, thus inter-valley backscattering can take 
place (if the potential is short range), since in this case the 
eigenvalue of h does not change sign. This discussion will be 
of importance for Sec. |V.B| 



in a magnetic field, the spinors change sign upon the 
transformation 9{k) — > 9{k) + 27r, as dictated by the 
spin-statistics theorem. 

The first strinking consequence of the chiral nature of 
electrons in graphene was the observation of the chiral 



quantum Hall effect ( Novoselov et al. 2005a Zhang et al. 



2005 ) , where the Hall conductivi ty is quantized as Oxy = 
2e^(l + 2n)/h, with n = 1,2, . . . (Gusynin and Sharapov 



2005 Peres et al. 2006 1 . The quantization rule follows 



Electron-electron interactions play no role in the half- 
integer or chiral quantum Hall effect. On the other hand, 
they are a crucial ingredient in the explanation of the 
fractional quantum Hall effect. During the first few years 
of graphene research, effects of electron-electron interac- 
tions have been elusive, but the recent observation of the 



1/3 fractional Hall plateau (Bolotin et al. 2009 Du et al. 



2009 Morpurgo 2009 ) , brings them to the forefront this 



active research area. It is a remarkable experimental fact 
that the fractional quantum Hall effect in graphene can 
be observed at magnetic fields of 2 T and persists up to 
a temperature of 20 K, for fields of 12 T. 

Using the results introduced above, we proceed to the 
discussion of several topics on electronic transport in 
graphene. 



III. CONDUCTIVITY AND CONDUCTANCE OF 
GRAPHENE AT THE DIRAC POINT 

As discussed in Sec. |TT] undoped graphene has its 
Fermi energy at the Dirac point, where the material has 
a vanishing density of states. This would naively sug- 
gest that the conductivity of undoped graphene should 
be zero. However, experiments challenge ones intuition 
and show a finite conductivity at zero energy (i. e., at 
the neutrality or Dirac point) . An example of a conduc- 
tivity curve of graphene is shown in Fig. |8] where we see 
that the experimental conductivity minimum, at V^ = 0, 
is of the order of ^ 4e^/ft. (horizontal dashed line). Val- 
ues of the conductivity minimum for several devices are 
given in Fig. [5] The existence of a conductivity mini- 
mum in graphene is also referred to as quantum-limited 
resistivity. 



from the nature of the Landau levels of Dirac electrons y\_ Sources of disorder 



(Johnson and Lippmannl 1949 Nieto and Taylor 1985 



Peres and Castro 2007 Rabi 1928 ) combined with the 



existence of the two valleys in graphene. 

The application of the chiral quantum Hall effect to 



metrology, in defining the resistance standard (Delahaye 



and Jeckelmann 20031, has clear advantages over the 



usual quantum Hall effect in the 2D electron gas, since 
the same experimental accuracy on the quantization of 
the Hall resistance can be achieved at higher tempera- 



As in any other metallic system, the electronic mobility 
in graphene is hindered by disorder. The sources of dis- 
order in graphene can vary, and can be due to adsorbed 
atoms (for example hydrogen) or molecules (for example 
hydrocarbons), extended defects, such as folded regions 
(wrinkles), vacancies, and topological defects [such as of 



Stone- Wales type, specially at the edges (Huang et al. 



2009 1]. Interestingly enough, in some particular cases, an 



tures (Giesbers et al. 2008 Poirier and Schopfe, 2010 extended defect in graphene can act as a ID conducting 



Tzalenchuk et al. 2010). At a temperature of 300 mK, channel (Lahiri et al. 20101. In addition, the system has 



the accuracy of the quantum Hall resistance quantiza- 
tion has been shown to be of 3 parts per billion, in 



a certain amount of rippling (random strain) ([Katsnelson 



and Geim 2008 Meyer et al.\ 20071, so it is not a perfect 



monolayer epitaxial graphene (Poirier and Schopfe 2010 planar lattice, and it has rough edges, which can exhibit 
Tzalenchuk et al. 2010). Also, the quantum Hall ef- scrolling (Fogler et al. 2010). We should note that, al- 
though the formation of vacancies is energetically unfa- 
vorable, the existence of adatoms and adsorbed hydrocar- 
bons is likely, originating from the isolation method and 
exposure to the environment. Such adsorbed atoms can 



feet in graphene has been observed at room tempera- 
ture (Novoselov et al. 2007[ ) and recently in epitaxial 
graphene as well (Wu et aL\ 2009[), which can be pro- 



duced in quasi- free standing form ( Riedl et al. 2009 1 



be imaged by transmission electron microscopy (Meyer 



et al. 2008). Additionally, the electrostatic random po- 



tential at the surface of the silicon-oxide substrate acts as 
an additional scattering source, originated from charged 



impurities (jZhang et al. 20091 



To a good practical approximation, an adsorbed hy- 
drocarbon, when binding covalently to the 2pz orbital of 
a given carbon atom of graphene, effectively removes the 
2pz electron from participating in the electric transport, 
by forming a cr— bond. Since the electron wave-function 
is spatially confined, the impurity can effectively act as 
a vacancy. This latter type of defects induce resonant 
states at, or close to, the Dirac point (see below). 

Another way of looking at this problem is to con- 
sider that, say, an hydrogen atom when binding cova- 



lently (Katoch et al. 2010 1 to a carbon atom in graphene 



changes locally the hybridization from pure sp2 to par- 
tially sp3 and creates, as before, a resonant impurity at 



that site ( Castro Neto and Guinea 2009 Robinson et al. 



2008). In this latter sense, both local potentials and 



adatoms have a similar effect (Stauber et al. 2008a I. The 



change of the chemical bonds from pure sp2 to partially 
sp3 adds an additional scattering effect originated from 



the enhancement of spin-orbit coupling ( Castro Neto and 



Guinea 2009) 



Combined with charged scatterers, the resonant scat- 
tering mechanism is currently ascending as one of the 
dominant processes limiting the electronic mobility in 



graphene ( |Ni et al.\ |2010| ) 



The resonant scattering mechanism is easy to under- 
stand by considering a simple model. We add to the 
Hamiltonian (111 a contribution from an impurity bind- 
ing covalently to a carbon atom at site i?„ — 0. Such a 
situation adds to the Hamiltonian a term of the form 
ifrs = (T^|ad)(A,0| + H.c.) + ead|ad)(ad|, where V is 
the hybridization between the adatom (or a carbon atom 
of a hydrocarbon molecule) and a given carbon atom of 
graphene, tad is the relative (to graphene's carbon atoms) 
on-site energy of the electron in the adatom, and |ad) 
is the ket representing the state of the electron in the 
adatom. Taking the wave function to be of the form |^) — 

E„[^(fl„)|Ai?n>+S(H„ + 52)|S,ii„+52)]+Cad|ad), 

the Schrodinger equation at the site i?„ — reads 

EA{Q) ~ yCad = -t[B{8i) + B{62) + B{8^)] , (5) 
{E - ead)Cad = VA{{)) . (6) 



Solving for Cad, we obtain 



-t[B[5^) + B{52) + B{5^)]^EA{Q) 



V^AjQ) 
E — ead 



(7) 



The resonant effect is included in the last term of Eq. 
(I?! , which represents a local potential of the form VeS — 
V^/{E — Ead)- Equation ^ contains two interesting 
regimes: (i) when \E\ <C Cad, the adatom acts as an ef- 
fective local potential of strength g^s = T^^/cad- If 5cff is 



large, the adatom acts roughly as an effective vacancy; a 
vacancy is characterized by goff = oo; (ii) when E sa e^d, 
the hopping from the carbon atom at position i?„ = 
to its nearest neighbors is suppressed [effectively we have 
t ^ {E — ead)i]i and the adatom acts roughly and again 
as an effective vacancy at energies close to ead- Therefore, 
either by inducing an effective local potential or by sup- 
pressing the nearby hopping we see that such mechanism 
increases the likelihood of an electron being trapped for a 
longer time in the vicinity of the adatom, thus generating 
a resonant state. 

If ead — 0, then the resonant states will be exactly 
at the Dirac point, and this is expected to happen for 
adsorbed hydrocarbon molecules. It is then the job of 
quantum chemical calculations to determine the value of 
the parameters ead and V (Robinson et al. |2008 , Wehlirig| 



et al. 2009 20101. Recently obtained typical values are 
1/ ~ 2i ~ 5 eV and ead ^ -0.2 ^ Wehling e^"al||2010[ ), lead- 
ing to gcff ^ 100 eV, a rather strong on-site potential. 
Finally, the calculation of the transport properties for 
such a model can be performed using the T— matrix ap- 



proach (Peres ei aZ. 2007a 2009b; Robinson eiaZ. 2008). 



Its derivation is elementary, using the simple model de- 
scribed above. It is well known that the T matrix for 



a local potential of intensity wo reads (Bena and Kivel- 

r 



son 



2005 Peres et al. 2006) T{E) = vo[l - vqG r{E)] 



Then, for an adatom we inust have 



T{E) 



Kff 



y2 



1 - V,sGr{E) ^ - ead - V^Gr{E) 



(8) 



Using Eq. (p|, it is simple to compute the transport 
relaxation time T{ep) at the Fermi energy ep using 
hlT{ep) = T:ni\T(ep)\'^ p{ep) , where n.^ is the concen- 
tration of impurities per unit cell. From the knowl- 
edge of T{ep), the conductivity of graphene follows 



Sec. 
ED 



from Boltzmann's transport equation (Ziman 1979) (see 
VI.C| . The function Gr{E) 

UW/D^)-i7T\E\/D^ 



reads: 
with D ~ 3t. 



Gr{E) 



It has been theoretically predicted that, in addition 
to their scattering effect, monovalent adatoms in diluted 
concentrations can create a gap in graphene's spectrum. 



by a mechanism called sublattice ordering (Cheianov 



et al. 2010 



Superlattices of vacancies (or adatoms) 



have the same effect (Martinazzo et al. 2010[) 



Midgap states ( [Jackiw 1984) are also produced by 



a model of pure vacancies (Pereira et al. 2006 2008), 



as shown in Fig. |4] and if a nearest neighbor hopping 
(t' ~0.4 eV) is included, the resonant states, while no 
longer at exactly the Dirac point, remain at energies close 



to it (Pereira et al. 2006 20081 



B. Calculation of the conductivity minimum for bulk graphene 
due to disorder 

It is certainly difficult to model all the different types 
of disorder just mentioned in a single calculation. We 
for the moment ignore this complexity and assume that 
electrons in graphene move in a random potential of the 
form V{r) = vq J2n=i H-^n ~ f) i where the position vec- 
tors Rn are random, vq is the strength of the potential, 
and Ni is the number of scattering centers. This model 
can be seen in the worst case scenario as zero order de- 
scription of the effect of impurities in graphene, although 



it has recently been used widely (Ando and Nakanishi 
T9981 lOstrovsky eTai] [20061 |2007| |Peres~f ai, 2006 



Suzuura and Ando 



2002 



Zheng and Ando 20021. In 



lU 



fact, for large vq this model mimics the resonant scat- 
terers physics. In what follows, we determine the conse- 
quences of the above random potential on the minimum 
conductivity of graphene. 

The usual approach to the calculation of the conduc- 
tivity uses the Kubo-Greenwood formula, obtained from 
linear response theory ( Mahan 2000 1 . The calculation 



proceeds in two steps ( Peres et al. 2006 Shon and Ando 



1998 Zheng and Ando 20021: first, the single parti 



cle Green's function in the presence of the disordered 
potential is computed in a self-consistent manner; sec- 
ond the current-current correlation function is obtained 
in terms of the single particle Green's function. This 
method is known as the self-consistent Born approxima- 
tion (SCBA). The final result of such calculation is a 
simple expression for the conductivity (j{£f) at the Fermi 
energy reading 



4e2 

TTIl 



(9) 



where K{tF) is a dimensionless function (Peres et al. 
20061 iStauber et aLl|2008a|); Eq. ([9| holds true both at 



finite Vq or when tig — >■ oo; from here on we consider this 
latter regime only. Since we describe the transport at the 
Dirac point, we need the value of K{ep) at zero chemical 
potential, which turns out to be K{Q) ~ 1. This result 
is essentially insensitive to the concentration of impuri- 
ties Ui = NiAc/A {A is the area of the sample and n^ 
is the concentration of impurities per unit cell) . The be- 
havior of K{0) as function of n^ is shown in the inset of 
Fig. |4] as stated, its value is 1. We have, therefore, ob- 
tained a universal value for the conductivity minimum of 
graphene Umin = 4e^/(7r/i) independent of the impurity 
concentration, even if the concentration of impurities is a 
small number. Many have reached the same result using 
different approaches (Dora et al. 2008 Ziegler 20071. A 



question naturally arises: How does one understand the 
result given by Eq. (|9l? To that end, we compute the 
density of states of disordered graphene. 

We have compared a calculation of the density of states 
as given by the SCBA with that given by an exact nu- 




Figure 4 (Color online) Numerically exact density of states 
(solid lines), in the limit i;o — >■ cx3, and the corresponding 
SCBA calculation (dashed lines), for the same impurity con- 
centrations rii. The maximum of the SCBA density of states, 
at S = 0, follows the rule p{E) ~ 0.2^/n7 eV~\ Inset: The 
function K{0) is plotted for a range of impurity concentra- 
tions spanning three orders of magnitude. (The numerically 
exact calculations are courtesy of Vitor M. Pereira.) 



merical method (Pereira et al 



2006 2008 



In Fig. |4] 



we show two sets of calculations for the density of states 
close to the Dirac point {E ~0). Two features are clear 
from these calculations. First, the disorder only affects 
the DOS close to the Dirac point, rendering it finite; sec- 
ond the SCBA introduces a smoothing of DOS around 
E ^0, but its value essentially agrees with the exact one, 
except at energies very close to E — 0. The finite den- 
sity of states close to the Dirac point is due to the wings 
of the resonant states forming at zero energy. The same 
behavior is seen in the local density of states around a 
single vacancy and in the corresponding scanning tun- 



nelling microscopy current (Peres et al. 2007a 2009b). 



The above comparison shows that the SCBA gives a 
reasonable description of the density of states close to 
the Dirac point, and this gives us a certain amount of 
confidence in the calculation of the conductivity cr{eF) 
based on the same approximation. A comment on the 
behavior of the numerical DOS close to zero energy is in 
order: the sharp feature at precisely E — seen in the 
exact numerical solution arises from the presence of zero- 
energy quasilocalized modes, induced by the vacancies in 
the lattice (IPereira et all |2006l 120081). These localized 



states are clearly not captured by the SCBA. 

In short, the finiteness of the conductivity at the Dirac 
point is a consequence of the finiteness of the DOS at 
E ~0 due to disorder, even when the concentration of im- 
purities is small, since K{0) is essentially constant over 
several orders of magnitude of impurity concentration. 
Furthermore, there is a strong criticism in the literature 
regarding the application of the SCBA approach to de- 



scribe the physics at ej? ~ ( Aleiner and Efetov 2006 1 , 



but not at finite ep, as long as weak localization effects 



are not important (see Sec. V.A|. 



C. Calculation of the conductivity minimum for pristine 
graphene ribbons 

Graphene ribbons have been produced by different 
methods: etching of exfoliated graphene (Han et al. 



20071, using chemical reactions (Jiao et al. 2010 Li 



E 
o 




et al] 2008a I, unzipping carbon nanotubes (Kosynkin 



et al. 2009 1 , and tailoring them by scanning tunneling 
microscope lithography (Tapaszto et al. 20081. Much 



of the experimental challenge regarding the production 
of nanoribbons is related to the discovery of an experi- 
mental procedure allowing, in a systematic way, the en- 
gineering of ribbons with fixed widths and perfect edges, 
together with a detailed characterization of their trans- 



port properties (Han et al. 20101. 



Figure 5 (Color online) The conductivity of graphene at 
the neutrality point. Left: conductivity minimum from the 
Manchester's group [Data from (Geim and Novoselov 2007 1]. 
Right: Conductivity minimum from the Columbia's group 
[data from ( Tan et al.\ |2007 l[ as function of the mobility 
of the devices. In both panels, several devices with differ- 
ent mobilities n were measured. The mobility is defined as 
fj, — a{eF)/{en), where n is the electron density, and is a mea- 
sure of the amount of disorder in the system. The constant go 

stands for twice the quantum of conductance, 2e^//i ~0.078 regime uses Landauer's formalism (jNazarov and Blanter 
(kn)"'^, and is represented by the dashed line. The dash- 
dotted line represents the value a = (Jmin, obtained in Sec. 



The previous section addressed the problem of the con- 
ductivity minimum of graphene from the point of view of 
disorder. Another relevant problem is that of the trans- 
port properties of pristine ribbons, where electrons can be 
in the ballistic regime. The problem we are about to dis- 
cuss is a rather interesting one, since electronic transport 
will proceed via evanescent modes, whereas in normal 
metals charge transport is associated with propagating 
states. 

We show below, and also in this case, the system has 
a finite conductivity, which in some conditions has the 
same value we found in Sec. |III.B| although the phys- 
ical mechanism is different. The approach to the cal- 
culation of the conductivity of ribbons in the ballistic 



20091, where the relevant quantity to be computed is 



the conductance of the system, which can formally take 
into account quantum interference effects, absent from 
the elementary Boltzmann's transport theory (but see 
Sec. 



V.Al 



If one considers that the scattering centers are charged 
impurities, the conductivity of graphene, at the neu- 
trality point, acquires the form (Fogler 20091 crmin — 
{e^ /h)cC with c = 0.5 ± 0.05 and C the solution of 
the transcendent equation C = \n[C/Aa'^^), where a'^^ 
is the effective fine structure constant of graphene (see 
Sec 



The measurements in bulk metals of dc-transport prop- 
erties allows one obtain directly the resistance R of the 
sample, from which the linear conductance G — 1/R can 
be determined. In bulk metals, we can define a mate- 
rial intrinsic quantity, the conductivity a. Taking the 
example of a 2D system, we have a = GL^/Ly, where 



IV.C| . This result for the conductivity minimum is L^ and Ly are the longitudinal and transverse dimen- 
sions of the bulk sample, respectively. The conductivity 
is a well-defined quantity whenever the system is large 
enough, such that the electronic current is homogeneous 
and insensitive to variations of the impurities' position 
from sample to sample. In this regime the transport is 
well described by Boltzmann's transport equation. The 
validity of this equation assumes that ( [Ferry et al.\\2009\ : 
(i) the scattering process is local in space and time, (ii) 
the scattering is weak and the electric field is small, and 
(iii) the de Broglie wavelength of the electron at the Fermi 
surface is much smaller than the distance between impu- 
rities. The systems amenable to such description are said 
to be self averaging. (In 2D, both a and G have the same 
units, 1/57.) 



different from that obtained for strong short-range scat- 
terers. 

Finally, the measured conductivity minimum, as shown 
in Fig. I5] has the same order of magnitude as that given 
by cTmin, but is larger than this value and has a finite 
variance. It is important to note that the assessment of 
the transport properties of graphene at the neutrality (or 
Dirac) point can be strongly affected by the used probe 
geometry, the use of invasive contacts, or the lack of ef- 



fective control on the sample's homogeneity (Blake et al. 
2009[ [Connolly et al. 2010 1. These effects are responsible 



for the differences in the two sets of measurements given 
in Fig. [5] 



When the system's size is reduced, we enter the realm 
of mesoscopic physics. It is instructive to compute the 
order of magnitude of the number of impurities in a 
graphene flake with an area of A = L? , and L = 0.25 
^m (see Sec. If]). Taking Ni/A = 5 x 10^^ cm^^ as a 
typical impurities' concentration in graphene (see Sec. 
IV. B for understanding the origin of this number), we 
obtain Ni ^ 3 x 10^ impurities. The typical distance 
between impurities is d ^ y/A/Nl ^0.02 /im; a typi- 
cal Fermi wave number for the electrons is graphene is 
(see Sec. 



(Katsnelson 2006 Tworzydlo et al. 20061. In the left 



IV I, from which follows that 



kp ~0.003 /im" 
de Broglie wavelength of the electrons at the Fermi sur- 
face is Xp = 2ti /kp ~0.02 /im, making d and Xp of the 
same order of magnitude. In this regime, the current be- 
comes non-homogeneous and sensitive to the position of 
the impurities in the material. Then, the conductance 
shows fluctuations from sample to sample, and the con- 
cept of conductivity loses its meaning. Metallic systems 
such as graphene are considered highly conducting but 
disordered metals. The behavior of the electrons becomes 
sensitive to the metal contacts, surfaces, and interfaces 
as well and quantum mechanical interference effects be- 
come important. Due to these interference effects, the 
transport properties of mesoscopic systems in the ballis- 
tic regime are better assessed by the Landauer's formal- 
ism (iFerry ei aZ.l|2009l. 



In calculating the conductance of pristine graphene rib- 
bons, we assume a ribbon of length L^ and width Ly , con- 
nected to heavily doped (say, with electrons) graphene 
leads (see Fig. 10 for the geometry of the device). The 



doped graphene leads will act as electron reservoirs, and 
the doping is modeled by gating the leads at a potential 

Since the leads are gated, there is a mismatch between 
the longitudinal momentum k^ of the electrons in the 
leads and in the central part of the device, where the 
undoped graphene ribbon lies; in the device electrons 
have longitudinal momentum q^ . The momentum ky is in 
this case a conserved quantity. The problem is then that 
of computing the transmission amplitude for an electron 
coming from the left lead to emerge at the right one. The 
energy of the electrons at the right and left leads is given 
hy E — —eVgztvp^/k'^ + k^', in the central region the en- 



ergy is given hy E ~ ztvp^/q^ + k"^ . We further impose 



periodic boundary conditions along the transverse direc- 
tion, which gives ky = iTmjLy with n = 0, ±1,±2, . . .. 
Since we are interested in graphene 's transport properties 
at the Dirac point, we have to consider the case of zero en- 
ergy. For this energy, the solution of \ q^ + fc^ = gives 

qx = iky and therefore the propagation of the electrons 
in the central region proceeds by means of evanescent 
waves. 

The scattering problem requires writing the wave func- 
tion on the left and right leads, and on the central region 



lead, the wave function, up to a multiplicative factor of 
g*yfe»^ reads 

^Lir) = (^ J^,,^ ) e^'===-+r„ (^ _^i,(,) ) e-"=^- . (10) 

In the central region the wave function can be written as 



ipcir) = a„ 






+ bn 



Finally, in the right lead we have 






(11) 



(12) 



The calculation of the transmission amounts to impos- 
ing the continuity of the wave function at a; = and 
X = Lx and determining the transmission amplitude i„ 
from which the transmission associated to a given trans- 
verse mode n is obtained as T„ = |i„p (to each quan- 
tized ky momentum corresponds a n transverse mode). 
The final result for the total transmission at zero en- 



E 



ergy is ( |Katsnelson[ |2006[ [Tworzydlo et al.\ |2006 [): T 



^„l/cosh (kyLx). As stated, the conduc- 



tance G is expressed in terms of the conductivity a as 
G — 4e^r/ft, — cLy/Lj.- In the regime Ly/Lx ^ 1, cor- 
responding to ballistic transport, we have T ~ Ly/{LxTr), 
and therefore a = a„iin, the same value obtained in Eq. 
(l9|, due to disorder. We stress that, for graphene rib- 
bons, only in the regime Ly/Lx 3> 1 is the conductivity 
a well defined quantity, since only in this case is this 
quantity independent off the aspect ratio of the ribbon. 
The extension of this type of calculations to finite tem- 
peratures is elementary, and it follows from the Lan- 
dauer's formalism as well. Such theoretical investiga- 
tions were done and the results seem to be in qualitative 
agreement with transport measurements made in high- 



20091 



mobility suspended graphene (Miiller et al. 

The conductance of ribbons, with aspect ratio 
Ly/Lx ^ 1, was experimentally measured, and the value 
c = Umin was obtained (Danneau et al. 



2008 Miao 



et al. . 2007 1 in agreement with the previous result. There 



are, however, difficulties associated with measuring the 
conductivity of graphene ribbons at the neutrality point 



(Blake et al. 20091, since inhomogeneous samples tend 



to overvalue the minimum of conductivity and two-probe 
measurements are generally expected to undervalue it 
(Blake et al. 2009 1 . Due to these subtleties, there is some 



reserve in the community ( Blake et al. 2009 1 regarding 



the measured conductances (Danneau et al. 2008 Miao 



et al. 20071. 



We note that the above result for CTmin, being equal to 
that computed in Sec. |III.B[ has a different physical ori- 
gin. The result obtained here is only valid in the regime 
Ly/Lx S> 1; when the system is in the ballistic regime. 



However, one must recognize that the presence of the 
evanescent modes in the above calculation produces a fi- 
nite density of states at the Dirac point, precisely what 
happens in the bulk disordered graphene calculation dis- 
cussed in Sec. |III.B[ When the calculation just described 
for graphene in ballistic regime (Ly 3> L^) includes the 
effect of resonant scatterers, the conductance is corrected 
by the value SG = Aamin/T^, per resonant scatterer (Titov 



(Titov 20071 



et al. 20101, that is, we have impurity-assisted tunnehng 



As a last comment, we note that the important top- 



ics of edge disorder (Gallagher et al. 2010 Lewenkopf 



et al. 2008 Mucciolo et al. 2009 1 and Coulomb block 



ade in graphene nanoribbons are not considered in this 
Colloquium, since they have been considered elsewhere 
(Dubois et al. 20091. A review on the effect of disorder 



on the electronic transport in graphene nanoribbons is 



also available (Mucciolo and Lewenkopf 20101. 



D. Puddles 

We now address the fact that the model developed in 
Sec. |III.C| for the conductivity of graphene at the Dirac 
point is somewhat simplistic, since it assumes the possi- 
bility to have graphene with exactly zero electronic den- 
sity a.t E — Q, the neutrality point. 

The physics close to the Dirac point is different from 
that at finite densities (to be discussed in Sec . |IV| as sug- 
gested by the data shown in Fig. [6j In Sec. |IV| we show 
that the electronic density n can in graphene on top of 
silicon oxide be externally controlled by a gate potential 
Vg and given by n ~ 7.2 x 10^° x Vg cm~^. According 
to this equation, the electron density can be tuned all 
the way down to zero by changing the gate potential. 
However, in Fig. |6] we show that the absolute value of 
the electron density never drops below its theoretically 
predicted value for Vg = 2 V. This experimental fact 
hints for different physics close to the Dirac point, where 
the system shows important charge-density fluctuations 
caused by the random electrostatic potential due to sub- 
surface charged impurities. We, however, stress that in 
suspended annealed graphene the electronic density can 
be made as low as ^ 10* cm~^, which corresponds to 
about a single electron present in a micron size device. 
Additionally, graphene's topography shows corrugations, 
which are probably due to roughness in the underlying 
Si02 surface and due to intrinsic ripples of the graphene 
sheet. 

As mentioned, the calculations in Sees. |III.B| and |III.C| 
assume as the starting point that graphene is a perfectly 
flat material, with null electronic density everywhere. 
However, experiments using a scanning single-electron 



transistor ( Martin et al 
models of Sees. HIT. Bl and 



2008) found that the idealized 
Ill.Cl do not hold. Those inves- 
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Figure 6 (Color online) Dependence of the absolute value of 
electron density on the gate voltage Vg. In Sec. |IV[ we show 
that the electron density n as function of the gate voltage 

Using this, the electron density 



for Vg = 3 V should be n : 
eV we should have n 



whereas for Vg = 2 
It is clear that for 



Vg = 3 the electronic density is the predicted one; however n 
never equals its predicted value for Vg < 2 V. Note that the 
vertical scale is logarithmic. (Data from M. Monteverde et al. 
(|Monteverde et al.^ 2010), courtesy of M. Monteverde.) 



to be a non-homogeneous system, with electron and hole 
puddles coexisting, with variations in the electronic den- 
sity in the range n G [—1,1] x 10^^ cm~^, which corre- 
sponds to a spatial variation of the surface electrostatic- 
potential in the range [—0.25, 0.25] V, with a full width at 
half maximum of 50 mV. The behavior shown in Fig. [Glis 
an indirect signature of this experimental fact. The ex- 
istence of puddles renders the descriptions of Sees 



III.B 



and nil. CI unsuitable. 

Posteriorly, a scanning tunneling microscopy (STM) 



study ( Zhang et al. 2009 1 was able to provide detailed 



information on the size and electronic density value of the 
puddles. This study allowed one to characterize the pud- 
dles with electron-density spatial resolution two orders 



of magnitude higher than previous investigations (Mar 
tin et al. 2008[). From the results of Sec. Ill] we can 



write a relation between the electronic density and the 
energy as n = E"^ / {TTVpfi^) . Given the presence of the 
puddles, the energy becomes function of position, as does 
the electronic density. We thus have a relation between 
the electronic density at the Dirac point and the energy, 
reading 



n{x,y) 






(13) 



The STM allows the determination of El,{x,y), from 
which n{x, y) is obtained. These studies revealed that the 
average lateral dimension of the puddles is of the order 
of ^20 nm [a theoretical study (Rossi and Das Sarma 



tigations ( Martin et al., |2008|) found undoped graphene 



|2008[ ) obtained a similar value], and that each of these 
puddles contains, on average, a charge of 0.3 ± 0.2 elec- 



10 



tron. A Kohn-Sham theory of the carrier-density dis- 
tribution of massless Dirac fermions in the presence of 
arbitrary external potentials has also predicted the exis- 



tence of the puddles (Polini et al. 20081. It was experi 



mentally determined that the topographical corrugations 
in graphene are about an order of magnitude smaller 
than the puddles' size, and therefore cannot justify their 
origin. Indeed, it was established that individual sub- 
surface charged impurities are responsible for the forma- 
tion of the puddles. It was estimated that the charge 
fluctuations associated with a single of those impurities 
is of the order of 0.07 ± 0.03 electron. There is, there- 
fore, a consensus that the physics of the puddles is due to 
charged scatterers. The origin of such charged scatterers 
is likely to be due to chemical species physisorbed onto 
graphene, which have been trapped in between the sub- 
strate and the graphene sheet during fabrication process 
of the device. 



From a theoretical point of view, graphene in the pud- 
dles regime can be thought as a random resistor net- 
work (Cheianov et al. 20071. Since Klein tunneling 



( Beenakker 2008 1 is exponentially suppressed if the bar- 



riers are not perfect potential steps ( Cheianov and Fal'ko 



20061, a large electronic transmission will not occur, ex- 



cept for perfectly normal incidence [things are markedly 
different for magnetic barriers as opposed to electrostatic 



ones (Martino et al. 2007)]; the essential physics relating 
the smoothness of potential barriers to the suppression 
of Klein tunneling was studied in the early days of rel- 
ativistic quantum mechanics, following a suggestion by 



Bohr (Christillin and d'Emilio 2007 Sauter 19311. The 



validity of the random resistor model depends on the 
assumption that transport is incoherent at scales larger 
than the puddle sizes. Due to Klein tunneling, massless 
Dirac electrons cannot also localize (Anderson localiza- 



tion) ( Cheianov and Fal'ko 2006 Katsnelson et al.\ 2006 



Lewenkopf et al. 2008 Mucciolo et al. 2009 ) under the 



effect of the random electrostatic potential (long-range 
scatterers) creating the puddles; this accounts for the 
finite conductivity of graphene at the Dirac point. As 
discussed in Sec. |V.B| long range scatterers preclude the 
possibility of weak localization effects, and since the elec- 
trostatic potential variations can be attributed to charged 
impurities, the description of transport at the Dirac point 
based on such type of scatterers seems to be the correct 



approach (Adam et all 20071. 



We note that intra-cone backscattering (see Sec. V.B ) 



has been shown to be present in graphene (Zhang et al. 



teresting experimental fact. Finally, when strong inter- 
valley scattering is present, electrons in graphene can lo- 
calize. 



IV. THE TRANSPORT PROPERTIES OF GRAPHENE AT 
FINITE ELECTRONIC DENSITY 

In the previous section we discussed STM experiments 



supporting the theory iPolini et al. 2008 Rossi and Das 



Sarma 2008 1 that charge scatterers dominate the elec- 



tronic transport of neutral graphene. In the ensuing sec- 
tions, we discuss transport in doped graphene, analyzing 
the role that resonant and charged scatterers play in this 
regime. 



A. The dependence of the conductivity on the gate voltage 

We now discuss the dependence of the conductivity 
of graphene on the gate voltage, considering two differ- 
ent types of scatterers: resonant scatterers (strong short- 
range scatterers) and charged impurities. We shall not 
discuss here scattering from random strain (Katsnelson 
and GeimI 2008), which we defer to Sec 



IV.D 



The electronic density in graphene can be controlled by 
the back-gate of a device engineered as a plane capacitor 
- a field effect transistor, made of silicon oxide (relative 
permittivity e = 3.9), with a thickness b of ^^300 nm. 
According to elementary electrostatics, the electric field 
in the dielectric is given by -Bcap = en/(eoe), with n the 
surface electronic density of graphene, which acts as one 
plate of the capacitor. The gate potential is related to 
the electric field by i^cap — Vg/b, and so the density of 
induced charge is n = eo^yg/i^b). Inserting the numerical 
values of e and b, we obtain n = aV„, with a ~ 7.2 x 10^'^ 



V~ -cm" . The Fermi momentum kp is obtained via 
kp — y/anVg, a result derived by counting the states in 
momentum space up to kp. 



Ever since the original paper on graphene (Novoselov 



et al. 2004), demonstrating the ambipolar field effect, it 



became clear that the conductivity of graphene depends 
on the gate voltage in some circumstances roughly as 
<j{ep) ex Vg] this is shown in Fig. [sl after some replotting 
of the data (solid curve on the right panel) . Experiments 
also show conductivities presenting a sub-linear behav- 
ior; see Figs, [t] (solid curves) and Fig. [s] (dashed curve 
on the right panel). Mobilities, a measure of the qual- 
ity of the electronic transport (see caption of Fig. [s] for 
the definition of the mobility /i), as high as /i ^ 1 x 10^ 
cm^-V~^-s~^, have been indirectly measured by Landau 



level spectroscopy ( Li and Andrei 2007 Li et al. 2009 1 



of graphene fiakes on top of graphite (Neugebauer et al. 



2009 1 , raising the question of how perfect can graphene 



be ( Neugebauer et al. 2009 1 . Ultimately, the answer re 



quires the identification of the limiting sources of elec- 
tronic scattering in graphene (among those listed in Sec. 



2009), which in view of Klein tunneling is a rather in- III. A I 



An approach combining Fermi's golden rule, Boltz- 
niann equation, the Coulomb potential created by 
screened charged impurities, and a random phase approx- 
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imation calculation of the dielectric function of graphene 
(Shung 1986 Wunsch et aL[p006 l gave a first good ac- 
count of the observed <j{eF) oc Vg behavior for graphene's 



conductivity (Adam et al.\ 2008 2007 Hwang and Das 



Sarma 2008). When graphene was doped with potassium 



(Chen et al. 2008), the measured conductivity agreed 
with the theory (Adam et al. 2007), as expected. [We 



note that the conductivity of graphene covered by metal 
clusters it is still far from being fully understood (Pi 



et al. 20091. 



Using the same approach for a delta-function potential 



(Adam et al. 2008 2007 Hwang and Das Sarma 2008), 



the prototype of a short range scatterer, the computed 
conductivity is a constant number, independent of the 
gate voltage and of the dielectric constant of the medium. 
In what follows, we argue that this result is inconsis- 
tent. We note that an attempt to solve the Lippmann- 
Schwinger equation for a delta-function potential showed 
that this problem is ill defined (regularization of the prob- 
lem is required in order to have a well-defined problem; 
as usual, this procedure introduces a length scale. This 
length scale is interpreted as the range of the short-range 
potential), and therefore the first Born approximation 
cannot be trusted. Indeed, exact numerical calculations 
show that the first Born approximation is inadequate 
to describe the role of strong short-range scatterers in 
graphene (Klos and Zozoulenko 2010). At the same 



time, a numerical calculation using the Kubo-Greenwood 
formalism ( Nomura and MacDonald") 2007[ ) showed that 
cr(ei?) oc Vg for charged impurities (the level broadening 
due to scattering was however introduced by hand) . The 
same work (Nomura and MacDonald 2007) also showed 



that short-range impurities do produce a conductivity 
that depends on the gate voltage, but in a sub-linear 
manner. Also, previous calculations of cy^ep) based on 
the SCBA showed that strong short-range scatterers, de- 
scribed by delta- function potentials, do give rise to a 



gate- volt age-dependent conductivity (Peres et al. . 2006 
Shon and Ando 1998), a result embodied in Eq. ([9|). A 



similar conclusion was obtained from a semi-classical ap- 
proach taking into account the chiral nature of massless 



Dirac fermions ( Trushin and Schliemann 2008 ) 



The two different results - those based on Fermi's 
golden rule, as opposed to those obtained from the 
SCBA, for strong short-range potentials - are easily un- 
derstood: the SCBA is a non-perturbative method, suit- 
able for strong short-range potentials, which takes into 
account the large deviation of the wave function, within 
the potential range, from the usual plane wave used in 
the first Born approximation, as pointed out by Peierls 



(Peierls 1979): indeed, the first Born approximation pro- 



duces a large scattering cross section, whereas the exact 
calculation gives a small value. Since the conductivity 
depends on the scattering (transport) cross section, an 
incorrect determination of it will give, at least, an incor- 
rect value for the impurity concentration in the material. 



Unfortunately, the reliance on the result based on Fermi's 
golden rule is widespread in the community and is being 



used to fit the experimental data (Hong et al. 2009), at 



the same time that the resonant scattering mechanism 
points toward the presence in the material of strong lo- 



calized potentials (see Sec III.A). 

If it is certain that some amount of charged impurities 
is present at the silicon oxide-graphene interface (respon- 
sible for the electron and holes puddles) , it is no less true 
that experiments do not rule out other sources of scat- 
tering. Indeed, recent experiments showed that both ad- 
sorbed hydrogen and vacancies led to conductivity curves 
indistinguishable in form from those of pristine graphene 
(iChen et all 120091 lEHas et all 120091). The presence of 



these short-range scatterers - vacancies and hydrogen - is 



signaled by a significant Raman D— band intensity (Chen 



et al. 2009), since they couple electron states from the 



K and K' valleys (see Sec. V.B I. By the same token, the 



presence in pristine graphene of such Z?— band would be 
the signature of the presence of short-range scatterers in 
the material. Detailed Raman investigations in pristine 



graphene have been carried out (Ni et al. 2010), showing 



that, indeed, a small D—peak is present in the Raman 
spectrum of the pristine material. Subsequent transport 



experiments (Ni et al. 2010) support strong short-range 



scatterers as the limiting source of scattering in graphene. 

An experiment especially designed to address the 

importance of charged impurities used devices with 



dielectrics having high permittivity constants (Pono- 



marenko et al. 2009 1 . These experiments did not ex- 



clude completely the contribution of this type of impu- 
rities, but did challenge the idea that charged impurities 
are the main source of scattering in graphene. 

On the other hand, in another set of experiments, an 
apparently similar investigation was done, but with ice 
layers on top of graphene and reaching a different conclu- 
sion. It was argued that the results were consistent with 
charge scattering dJang et al. 20081. There is, however. 



at least one difficulty with the arguments developed in 
that work: the number of ice atomic layers was at the 
most six and therefore can hardly be considered an infi- 
nite dielectric made of ice; the lines of the electric field 



are essentially in the vacuum (Silvester and Ferrari 1996 



Sometani 2000). 



A number of questions can still be asked ( Monteverde 



20071: 



et al. 2010 Ponomarenko et al. 2009 Schedin et al. 



1. In a study of graphene's sensitivity to gases 



( [Schedin et aZ.||2007p (NO2, H2O, and iodine acting 
as acceptors, whereas NH3, CO, and ethanol act- 
ing as donors), chemically-induced charge-carriers 
concentrations as large as 50 x 10^" cm~^ were 
achieved. The induced chemical doping shifted 
only the neutrality point of the conductivity curves, 
without any significant changes either in the shape 



12 



of those curves or in the mobility of the devices; the 
estimated concentrations of added c harged scatter- 
ing centers was high as 10^^ cm^^ (Schedin et al. 



and sub-Hnear behaviors can be accommodated within 
a model based on what is now called resonant scatter- 
ers dominated conductivity, giving rise to mid-gap states 



2007 1 . Why is it that no appreciable changes in ( Basko 2008 Robinson et al. 2008 Stauber et al. 2007 



the mobility were measured in these experiments? 
[One possible way out can be envisioned: the chem- 
ical dopants may cluster, and this would reduce the 



effectiveness of their scattering effect (Katsnelson 



et al. 2009 1.1 We also note that the definition of the 



mobility used in the analysis of the data ( Schedin 



et al. 2007) has been criticized in the literature 



(Hwang and Das Sarma 20071 



2. In a study designed to test the prediction (Hwang 



and Das Sarma 2008 1 of the charge scattering 



model for the ratio of the transport scattering time 
r and the elastic scattering time t^. for both mono- 
layer and bilayer graphene, the experiments found 
disagreement between the predicted behavior and 
the measured data, for both graphene systems. The 
measured deviations were found to be stronger for 
bilayer graphene (Monteverde et al. 20101. Fur- 
ther, it was found that the measured data agree 
with the resonant scattering mechanism. How to 
reconcile this set of measurements with models ex- 
plaining the mobility of both monolayer and bilayer 
graphene based on the charge scattering mechanism 
dZhu et al.\\2009\ 7 [We mention that Monteverde 's 
et al. results must be confronted with those of a 



similar experiment (Hong et al. 2009 ), reaching dif- 
ferent conclusions.] 

3. Since screening is strongly dependent on the value 
of the permitivity e of the surrounding medium, 
why is the mobility almost insensitive to changes of 
this parameter? For example, e for ethanol changes 
from 25 to 55 as the temperature drops from 300 
K down to ^160 K, but an experiment done in 
ethanol showed no variation of graphene 's mobility. 
We, however, note that in some experiments ( ^Pono-| 
marenko et al. 2009 1 a certain amount of variation 



in the mobility was measured in some devices upon 
changing the dielectric constant. This result does 
show that charged impurities play some role as scat- 
tering centers but apparently not the limiting one. 

The answers to the above questions remain debatable 
to some extent. The clarification of some of these issues 
could be taken to an ultimate test using a solid dielectric 
with a high relative permittivity. It just happens that 
strontium titanate (SrTiOs) has a relative permittivity 
of about 10,000 below T = 50 K, which suddenly drops 
to 300 when the temperature rises above 50 K. A device 
using such a dielectric would produce a dramatic change 
of the mobility upon a drop in temperature below 50 K. 

If we now refocus our attention on the role of strong 
short range scatterers, we recall that both the linear 



Wehling et al. 2009 20101, plus the additional effect of 



charged impurities, which, however, do not play the cen- 
tral role. On the other hand, the simplest model based 
on short range scatterers, in which the effect of charged 
impurities is ignored, does not account for the observed 
dependence of the mobility on the dielectric constant of 
the device, which has been shown experimentally to be 
present to some extent (Ponomarenko et al. 20091. In 
Sec. IV. B we present the main results of such a simple 



model, and in Sec. |IV.C| we include the role of charged 
impurities, and an improved model taking into account 
both types of scatterers is given. This latter model is 
a simple combination of results already available in the 
literature, albeit presented with a different emphasis. 

We stress that from the analysis of the SCBA results 
we see that features showing at energies close to the Dirac 
point are all proportional to y^ (recall that n^ is the 
density of impurities per unit cell), as shown in Fig. 
m This introduces an energy scale emin ^ fiVF^^/o-o, 
below which electron scattering based on plane waves 
breaks down, meaning that close to the Dirac point the 
results of Sees. |IV.B| and |IV.C| are expected not to hold. 

Both models based on mid-gap states or on charged 
impurities, presented below, fail to give a satisfactory 
account of the physics close to the Dirac point, since they 
are based on the scattering of plane waves. 

Finally, in suspended graphene (Bolotin et al. 2008a|b 



Du et al. 2008), where the material is hanging over 



a trench, mid-gap states are expected to survive, since 
some fraction of the corresponding scatterers will still be 
present, whereas charged impurities are expected to be 
absent. In the suspended situation, the mobility of the 
material will be limited only by resonant scatterers plus 
ripples induced by strain due to the electric field created 



by the gate (Bao et al. 2009 Fogler et al. 2008) 



B. Partial-wave description of resonant scatterers 

We assume the presence of short-range scatterers, 
which we model here as disks of radius R, and whose 



origin was discussed in Sec. III.B The effect of the scat- 
terers is such that the electron wave function is zero for 
r < R. The sizes of these disks are of the order of the size 
of the primitive lattice vectors. The circular shape takes 
the isotropy of the scattering process into a account, and 
the boundary condition allows a simple analytical solu- 
tion. A vacancy is one of the possible physical realiza- 
tions of the model just introduced. 

The Dirac Hamiltoniandslin polar coordinates r and 
ip reads (jHentschel and Guinea 2007, Peres et al. 2009a 
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Recher et al. 2007 1 



Hk = —ivph 



L_ 

L+ 



(14) 



with L± = e^'''^{d/dr±ir ^d/d(p). A particular solution 
of Eq. ( 14 1 with eigenvalue vphk has the form iprnifk) oc 

(J™(fcr)e-™'^,-^J™+l(fcr)e"'^'"+'^'^)^ where J^{z) is 
the regular Bessel function of first kind and integer order 
m. This solution corresponds to a partial wave in the 
angular momentum representation of the plane wave. In 
the presence of the potential created by the disk we write 
the trial wave function as 



1/ h\- Af J.n{kr)e^"'v 



-B 



iy„+i(A:r)e'(™+i)'^ 



(15) 



where Yra{z) is the irregular Bessel function of first kind 
and integer order m. We consider that at r = i? the 
wave function satisfies the zig-zag boundary conditions 



(Akhmerov and Beenakker 2008 Dominguez-Adame 



1990), i>i,m{kR)/il^2,TnikR) = 0, where ipi^mikr), 



with 

i = 1,2, is the j— component of the Dirac spinor. This 
boundary condition makes sense since, as discussed in 
Sec. |III.A[ a resonant scatterer can effectively behave 
as a vacancy; in turn, a vacancy is a three-site zig-zag 
edge. This boundary condition also represents the lim- 
iting case where the electronic probability fiux is zero 
through the region where the potential is finite. The 
phase shift 5m{kR) of the m partial wave is given by 
Hentschel and Guinea 



( Basko 



2008 



2007 



Katsnelson 



and Novoselov(|2007||Stauber et aLH2067||WehUng et al. 
2009| 



tan(5„j(fci?) 



Jmykii) 
YUkR) ' 



(16) 



The relative importance of the several phase-shifts to the 
transport scattering cross section axikR) depends on the 
value of kpR. In 2D. the differential cross section, 17(1^9), 
reads a{ip) = \f{(p)\'^, with f{if) given by 



rn— — 00 



^iipm^iS^ik) ^ 



sin[<5„(fc)] . (17) 



When kpR < 1, the s— wave phase shift SQ{kpR) is 
the dominant contribution. Making use of the relation 
l/r(fc) = n^VF(7TikR)/Ac ( |Ziman[ [T979| , where T(fc) is 
the transport relaxation time (see also Sec. V.A), and 
since the total transport cross section arikR) is obtained 
from 



arikR) 



27r 



d(pa{(p){l — COS (fi) , 



(18) 



the conductivity crikp) is given by (both spin and valley 
degeneracies included) 



(j{kp) — e V 



2 2 ^(fc) f N 2 P(eF) 

Ac niaT[kFR) 



(19) 



with (TTikpR) = 4sin So{kpR)/kF- Since we assume 
kpR < 1, we also have l/SoikpR) — 2hi{kFR)/'iT, and 
the final result for the dc-conductivity is then 



aikp) ^ 90-. ^In {^aTrVgR). 

An Ui 



(20) 



The result ( 20 1 for the conductivity holds as long as the 
Fermi momentum is larger than kp ^ emin/(^f^F) (re- 
call previous discussion) . The conclusion is that resonant 
scatterers (strong short-range impurities), giving rise to 
mid-gap states, give a conductivity that is gate voltage 
dependent, with sub-linear or quasi-linear dependence on 
Vg, depending on the size of the scattering disk R. Fur- 



thermore, the conductivity (20 1 is not independent of y^ 



even if we take R to be of the order of the carbon-carbon 
distance, ag. 

In the first experimental study of the conductivity of 



suspended graphene (Du et al. 2008), it was shown that 



this quantity is well described by fitting it to a model of 
mid-gap scattering states, that is Eq. (20). 



In Fig. |7] we fit the conductivity data of suspended 
and non-suspended graphene using Eq. ( 20 1 . In all cases 



a good fit is obtained, using impurity densities ranging 
from 1.3 X 10^^ cm~^ for non-suspended samples down 
to 0.7 X 10^" cm~^ for suspended ones. We note that 
for the devices termed KYI and K12 in Fig. [7] the rii 
values used are n^ ~ 0.8 x 10^^ cm~^ and n, ~ 1.3 x 10^^ 
cm~^, respectively. The attempt to fit the same data 



with charged scatterer s (Adam et al. 2008 2007 Hwang 
and Das Sarma 20081 gave concentrations of 2.2 x 10^^ 
cm~^ and 4.0 x 10^^ cm~^ for Kll and ^^12, respectively. 
The two set of numbers for Ui have the same order of 
magnitude, but Eq. (20) gives a better fit to the data. 



except for the measurements shown in the top left panel 
of Fig. [7] For this particular device, made of suspended 
graphene, the maximum of the measured conductivity is 
about two times smaller than that measured in the other 
suspended device, whose data are shown in the top right 
panel of Fig. [7] suggesting that additional sources of 
disorder may have been introduced during the fabrication 
process. We stress that the fits done in Fig. [Tj use the 
density of impurities as the only fitting parameter since 
R has to be of the order of Og ( Wehling et al.\ |2010 l; 
changes of order one in the value of R give small changes 
for the impurity concentration. 



Some phenomenological approaches (Morozov et al. 



2008 1 have tried to reconcile the measured sub-linear 



behavior of the conductivity of graphene with the lin- 
ear behavior (upon Vg) predicted by the charge scatter- 
ers model. In the case of the data given in the right 



14 




-80-60-40-20 20 40 60 80 -60 -40 -20 20 40 60 

gate voltage (V) gate voltage (V) 



Figure 7 (Color online) The conductivity of suspended and 
non-suspended graphene as function of the gate voltage. Top 
panels: Conductivity of suspended graphene from two dif- 
ferent devices, after current annealing {£ ~1 fiin). The 
top left panel corresponds to a device with fj, ~ 170, 000 
cm'^-V^^-s"'^, at electronic density of n : 



2x10" cm"^ [data 



from K. I. Bolotin et al. (Bolotin et al. 2008aj, courtesy of 
K. I. Bolotin]. The top right panel corresponds to a de- 



200, 000 cm' 

„„-2 



•V" 



^, at electronic d ensity of 
( Bolotin 



vice with fi 

n — 2 X 10^^ cm~^ [data from K. I. Bolotin et al 
\et al.\ |2008b| ), courtesy of K. I. Bolotin[. Bottom panels 
Conductivity of graphene on top of silicon oxide, correspond- 
ing to devices with fj, ^ 10,000 cm^-V'^^-s"^ (Data from Y.- 
W. Tan et al. {Tan et al\^007\ , courtesy of P. Kim) In all 
panels the fits use the model developed in Sec. (IV.Bl, with 
rii as the only fitting parameter; R has to be of the order of 



ao (Wehling et al 20101. In the legends, the concentration 



of impurities is per unit cell. 



panel of Fig. [7] a sub-linear behavior of the conduc- 
tivity upon Vg (for graphene on top of Si02) is evi- 
dent. A linear behavior could be recovered by defin- 
ing the measured resistivity ^measured as a sum of two 
terms ^measured = Pg + Ps, whcrc ps is a fitting param- 
eter. The conductivity aguh = 1/Pmoasured is sub-linear 
in Vg, whereas the conductivity crii„ = l/(pmoasurcd - Ps) 
shows linear behavior. The fitting parameter ps was as- 
sumed to be independent of Vg, and was attributed to 
short-range scatterers. The discussion presented above 
for strong short range scatterers showed that this type of 
disorder does produce sub-linear behavior of the conduc- 
tivity upon Vg, even in the case R ^ ao, with the same or- 
der of magnitude for impurity concentration as those pro- 



posed by the charge scattering mechanism (Adam et al. 



2008 2007 Hwang and Das Sarma 20081. It is obvious 



then that, using the curves calculated with the mid-gap 
state mechanism, it is still possible to obtain a linear de- 
pendence of the conductivity on Vg by assuming a ps 
fitting parameter as done in the phenomenological ap- 
proach (Morozov et al. 20081. From the discussion in 



this section it is fair to say that the origin of ps still 
needs clarification. Moreover, we can even ask the ques- 
tion whether the parameter ps is really needed for the 
interpretation of the data. 

Finally, we note that in our calculation we have not in- 
cluded the effect of inter- valley scattering, which is known 
to be present when the scatterers are short-range. If that 
effect is included, a contribution to the conductivity of 



the form given by Eq. ( 20 ) is found, albeit with a differ- 



ent numerical prefactor (Ostrovsky et al. 20061. 



C. Partial-wave description of Coulomb scatterers 

We now derive the contribution to the conductivity 
of graphene due to Coulomb scatterers. We can think 
of three alternative scenarios for the origin of Coulomb 
scatterers: either they exist independently of the reso- 
nant scatterers, or the latter can themselves be charged, 
carrying a fraction of the unit charge, or both cases can 
coexist. Experiments aiming at studying in detail the 
Raman D peak of pristine graphene, can also shed light 
on this aspect, by studying suspended pristine graphene 
before and after annealing. 

The solution of the Dirac equation in 2D for the 
Coulomb potential was obtained more than ten years 



ago (Lin 19991 for the sub-critical regime (see below), 
and rediscovered in the context of graphene by different 
groups (Novikov 2007a|b Pereira et al. 2007, Shytov 
et al. 2007 1, who also solved the case of the super-critical 



regime. 



It was known for some time (Lin 19971 that the 



Coulomb problem in the 2D Schrodinger equation pro- 
vides a total cross section which does not coincide 
with the first Born approximation. The same happens 



with the 2D Dirac equation (Lin 



Pereira et al. 2007 Shytov et al 



1999 



Novikov 2007a 



20071. The discussion 



given above, restrains us from accepting results based on 
the first Born approximation without a critical analysis. 
For the Coulomb problem, the Hamiltonian has the 
form H — HK + IZe^/{4:'KeQr) ( J a 2 x 2 identity matrix) 
and the solution is sought in the form 



V'j(r) 



/j(rfc)e''^(j'^i/2) 
±igj(rfc)e*'^(j+i/2) 



(21) 



with j a half-integer number. When the trial wave func- 



tion (21 1 is inserted into the Dirac equation we get 



e + g/r -dr - j/r 
dr - j/r e + g/r 



±igj 



(22) 



where j = m — 1/2, e 

2 



Zag, and ag 



e" /{A-Ke^vph) 



E/{vph), E is the energy, g — 
2.2 is graphene's fine 
structure constant. The solution to Eq. (22) has been 
obtained by several, and the central quantity is the phase 
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shift of the j partial wave, as in the case of Sec. 
For this problem the phase shifts read 

r{s + l + tg) 

with the property 5j{g) = 5-j{g), s = ^f - g'^ 
T{x) the usual gamma function; the sub-critical regime is 
defined by the condition g < 1/2. Contrary to the short- 
range scatterer problem, solved in Sec. |IV.B[ the phase 



(23) 



and 



shifts ( 23 ) do not depend on the energy of the incoming 



particle, but they do depend on the sign of the Coulomb 
potential and on which type of particle, an electron or 
a hole, is being scattered. The independence of 5j{g) 
on the energy is simple to understand from a straight- 
forward argument based on dimensional analysis: on the 
one hand, the Coulomb potential has no intrinsic length 
scale, and on the other hand, the particle's mass is null. 
These two facts show that the problem as a whole has 
no intrinsic length scale (some sort of Bohr's radius, as 
in the non-relativistic theory of the hydrogen atom) and 
therefore dimensionless numbers involving the momen- 
tum k cannot be formed, leading to the conclusion that 
the dimension of the cross section (dimension of length, 
in 2D) can only come from the momentum itself. That is, 
we are then bound to have (7T{k) oc 1/fc, which gives the 
linear dependence on the gate voltage. The electron-hole 
asymmetry of the cross section can partially account for 
the measured asymmetry of the conductivity curves. An- 
other source of electron-hole asymmetry of the conduc- 
tivity is originated in the metal contacts of the transistor 
( [Huard et a/., 2008 J . 

The above solution assumes that graphene is floating in 
vacuum. In a real experiment, graphene is on top of a di- 
electric, Si02 being the most common. Other dielectrics 
have also been used ( Ponomarenko et al. 2009 1 . In these 



experimental conditions, the value of g is different from 
that given above. 

In the case Z = 1 and for graphene on top of a di- 
electric, a charge e in between the dielectric (with rel- 
ative permittivity e^) and graphene, behaves effectively 



as a charge with a value ( Landau et al. 1984[ Slater and 



Frank 19691 of e^ = 2e/{l + e^). Additionally, the rel- 



ative pertimittivity of graphene due to electron-electron 
interactions is renormalized to er = 1 -|- 27rag/[2(l -I- ed)] 



( |Adam et am2007| [Gonzalez et ami999| |Shung||1986| ) 



These two effects combined give an effective fine struc- 
ture constant for graphene of 



a 



off 



Cr 1 + Ed 



(24) 



Using the same procedure of Sec. |IV.B[ the conductivity 
of graphene due to Coulomb scatterers reads 



e'^vp 



kEPJe-p) 
2niA{g) 



.90 



Sn^A{g) ^ 



ial 



(25) 



Dielectric 



Cd a„ 



A(g) A(-ff) 



H2O 



80 0.05 0.013 0.012 



Ethanol (160 K) 55 0.07 0.027 0.022 

Ethanol (300 K) 25 0.13 0.10 0.07 

Hf02 25 0.13 0.10 0.07 

Si02 I 4 0.37 1.20 0.46 

Table I Dependence of A{g) on the type of dielectric for both 
electrons, g > 0, and holes, g < 0. The impurities are assumed 
to have valence Z = —e. For ethanol, the dielectric function 
depends on temperature, as indicated between braces. Since 
graphene is in the ultra-relativistic limit, the value for A{g) 
cannot be obtained from adding only few partial waves. 



where A{g) = J2Z=~oc sin^('^m+i/2-<5m-i/2)- H is worth 
stressing that A{g) is different for particles, g > 0, and 
holes, (7 < 0, a behavior not captured by the first Born 
approximation (Adam et al. 2007). In Table IT] we give 
the numerical values for the quantity A (5), considering 
graphene on top of or submerged in different dielectrics. 
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Figure 8 (Color online) Resistivity and conductivity of 
graphene on top of silicon oxide. Top left panel: Raw data 
of a measurement of the resistivity, pmoasurod, of an exfoliated 
graphene sheet [data from S. V. Morozov at al. (|Morozov| 



et al.\ 2008)]. Right panel: Fit of the conductivity using 
Eq. (|20[ ) for the case where anub = 1/pmcasurcd, and using 
Eq. ( |25[ l for the case where aun — l/(pmcasurcd — ps), with 
PS = 100 Q. [th e value of ps = 100 Q is that used by S. V. 
Morozov at al. (Morozov et al. |2008[ )[. 

Data CJ"lin = t/ (Pmeasurcd 



Bottom left panel: 

ps) fitted with the conductivity 
formula given by Eq. (25 I (squares). For comparison, we give 



the theoretical conductivity curve (triangles) considering that 
the experiment had been done using Hf02 as a dielectric. This 
allows one to compare the modification of the numerical val- 
ues of a due to a substrate change. (Data from S. V. Morozov 
at al. ([Morozov et a/., ,2008 j, courtesy of A. K. Geim.) 
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We have now developed all the tools needed to perform 
the analysis of the data of Fig. [8] We start by fitting the 
data using the two models presented above separately, 
that is Eqs. (20) and (251. We must stress that each of 
these two models have only one fitting parameter: the 
concentration of impurities (in the model for resonant 
scattering, the parameter R is fixed by the size of the 
primitive cell). In the right panel of Fig. Islwe plot the 

data CTsub = 1/Pmoasured ^ud (Tiin = l/(Ancasurcd " Ps) 

using the raw data Pmcasurcd (given in the left top panel 
of the same figure) . In the case of CTgub , the data can be 
fitted using Eq. ( pO| ) for mid-gap states. A perfect fit 
to both negative and positive gate voltages is not pos- 
sible, since by construction the model developed in Sec. 
I IV. B I preserves electron-hole symmetry; an improvement 
which does not preserve electron-hole symmetry is easy 
to develop by considering a large finite value (as opposed 
to an infinite value) for the effective potential ^eff as dis- 



cussed in Sec. III.A| (see also (Araujo and Peres 2006 



Stauber et al.\ 2008a|); we show below that charged scat- 



terers can account for the loss of electron-hole symmetry 
of the conductivity curves as well. 

In the case of the data computed as crnn we fit the pos- 
itive gate voltage region with Eq. ( 25 1 for charged scat- 
terers. Note that since K{g) ^ A(— g) the computed con- 
ductivity has no electron-hole symmetry, an effect seen 
in the experiments ( Chen et al.\ 2008[ ). Nevertheless, al- 
though Eq. ( 25 I does break electron-hole symmetry, the 



magnitude of the computed effect is far too strong, and 
therefore Eq. ( 25 1 is not able to fit the data over the 

by assuming a single 
It is worth noting 



negative and positive range of Vg, 
concentration of charged scatterers. 
that the concentration of impurities used to fit a^uh and 
tJiin is essentially the same for both types of scatterers. 
In the left bottom panel of Fig. |8] we depict the con- 
ductivity values (triangles), had we performed the same 
experiment using Hf02 as a dielectric. To understand 
such a large change, we look at Table |T] where we show 
that K{g) can be reduced by one order of magnitude (pos- 
itive g) from Si02 to Hf02, leading to the large increase 
in the conductivity shown in Fig. [8] 

We now take into account, in a single model, the ef- 
fect of both strong short-range and charged scatterers. 
Computing the conductivity as Usub — 1/p, such that 
(Matthiessen's rule) 



1 



1 



P = 



f^shoit *^Coulomb 



(26) 



and with (Tghort computed using Eq. ( pO| and crcouiomb 
determined from Eq. (25), we can fit the data of Fig. 
[8] quite accurately, as shown in Fig. [9j In this figure, 
the concentration of impurities leading to mid-gap states 
used in the fit was rishort = 1-5 x 10^^ cm~^ and that 
for charged scatterers was ricouiomb = 2.4 x 10^" cm~^, 
about seven times smaller than nshort (these two concen- 
trations are, essentially, the only two fitting parameters 



in the model). The combined contributions from reso- 
nant scatterers and charged impurities allow one to fit 
CTsub over the whole Vg range, using a single value for the 
impurity concentrations. 
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Figure 9 (Color online) The conductivity data (Jsuh = 
l//9mcasurcd, also plotted lu Fig. p] but fitted using Eq. (26 1, 
which combines the effect of resonant and charged scatterers. 
In the legend, righort and ncoioumb refer to the concentration 
of resonant and charged scatterers, respectively. The concen- 
tration of impurities is the only fitting parameter used in the 
theory. [Data from S. V. Morozov at al. (Morozov et al 
2008[), courtesy of A. K. Geim.)] 



In conclusion, we made a thorough analysis of the role 
of resonant scatterers (which give rise to mid-gap states) 
and charge scatterers in the conductivity of graphene, 
and showed that a coherent picture emerges from a scat- 
tering analysis of the transport based on the exact cal- 
culation of the phase shifts of scattered chiral Dirac 
fermions, as opposed to a calculation based on the first 
Born approximation. 

Finally, we note that fine tuning details coming from 
the dependence of the dielectric constant of graphene on 
the wave vector (the polarization contributions) were not 
included in our simple model, except for the important ef- 
fect of the renormalization of the fine structure constant, 
which corresponds to the large wave number limit. 

Experiments will decide which scenario regarding the 
limiting source of scattering in graphene actually pre- 
vails. 



D. Transport across a strained region: A way of generating a 
transport gap 

The main limiting factor of all graphene properties, 
in what concerns its application to nanoelectronics, is, 
most likely, the lack of a true band gap, as opposed to 



the biased graphene bilayer ( Castro et al. 2007 McCann 



2006 McCann and Falko 2006). This fact can, however 



be overcome by creating a transport gap. 

In nowadays nanotechnology, understanding the effect 
of strain on the properties of devices is an essential step 
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toward the improvement of their performance. For exam- 
ple, characterizing how strain can improve the properties 
of silicon-based devices is a mainstream research topic 



(Hytch ei aZ. 20081. As stated in Sec. III. A both rip- 



ples and wrinkles can act as scattering centers as they 
effectively create random strain in the material, leading 
to a modification of the hopping energy t. In what fol- 
lows, we show that strain in graphene gives rise to a rich 
structure in the electronic and transport properties of the 
material. 



Being a 2D flexible membrane (Booth et al. 2008 Kim 



and Castro Neto 2008), stretching (which in graphene 
can be as large as 20%, being reversible) and bend- 
ing graphene in a controlled way is feasible (jFerralis 
|Kim e^ aZ.[ |2009| [Mohiuddin et al] \2009\ 



et al. 2008 



with consequences to the electronic (Pereira et al. 2009) 



and transport properties of the material (Fogler et al. 



2008; G uinea et gEj [20Tol [Pereira and Castro Neto| 12009 



Teague et al. 2009). As we see below, strain can be 



modeled by a fictitious gauge field (de Juan et al. 2007), 



which can then act as an effective magnetic field. In some 
circumstances, it was predicted that this effective mag- 



netic field can have an intensity as high as 10 T (Guinea 



et al. 2010 1, leading to a pseudo-magnetic quantum Hall 



effect. The presence of such an odd quantum Hall ef- 
fect can, in principle, be experimentally observed using 
scanning tunneling microscopy, which is a direct mea- 
sure of the density of states, and, therefore, sensitive to 
the reorganization of the spectrum due to the presence 
of the gauge field. Such type of experiments have been 



performed (Levy et al. 2010) and found strain-induced 



pseudo-magnetic fields greater than 300 T. 



In the case of suspended graphene (Teague et al. 



2009) there are two sources of strain. One is induced 



by the electric field produced by the gate, which pulls 
the graphene membrane downwards. The solution of 
the corresponding elasticity problem produces an effec- 
tive model where the effective vector potential is constant 



(Fogler et al. 2008 Pereira and Castro Neto 2009), pre 



cisely the model we discuss below. The other source of 
strain depends on the thermal properties of graphene. 
Graphene's thermal expansion coefficient is anomalously 



large and negative (Bao et al. 2009 Faugeras et al. 



2010), a feature which can be exploited to induce ID 



and 2D ripples (with a periodicity of about 300 nm) pos- 
sibly leading to novel strain-based engineered graphene 
devices ( Guinea et al.\ 2010 Pereira and Castro Neto 
2009|. 



Consider a graphene-based device where the central 
part of the material, of length L, is a graphene rib- 
bon under strain, with armchair edges oriented along the 
x— axis, as shown in Fig. 10 The strained part is then 



connected to two pristine leads. It is well understood that 
the effect of strain can be included in the Dirac Hamilto- 



nian in the form of a fictitious gauge field (^Castro Neto 



et al. 2009 Guinea ei aZ. 2010). The emergence of the 




Figure 10 (Color online) Scheme of a device made of strained 
graphene. The central region is the strained part. In the as- 
sumed model, the effect of strain is to modify the hopping 
connecting a given carbon atom to its three neighbors in such 
a way that two of the hoppings are equal. Shown is the mo- 
mentum of the electrons in the leads and in the strained re- 
gion. 



fictitious gauge field is simple to understand. We assume 
that the hopping along the armchair edge is modified rel- 
atively to its pristine value i as i — > i + Ai [a detailed 
study of how strain changes the value of the hopping was 
done using ab initio methods (Ribeiro et al. 2009j)]. This 



adds a term to the Hamiltonian, Eq. M, of the form 



At}^i\A,Rn){R„ 



aoU3:,B\ + H.c. 



(27) 



Passing from the tight-binding description to the contin 
uous model, the contribution from Eq, 
value at a given point in space 



(|27|) has a finite 
couples 



Since Eq. (|27| 
the sub-lattices A and B at the same point in space its 
contribution to the effective Hamiltonian, Eq. (l3|, has 



the simple form ( Castro Neto et al. 2009 1 



A{x, y)cTy = -e{x)e{L - x) UyUy , (28) 

evp 

implying that the Dirac Hamiltonian maintains its orig- 
inal form, but with p replaced p — > {px,Py + eAy). 
Clearly, we have a new Hamiltonian where the electrons 
now couple to a fictitious vector potential A through the 
usual minimal coupling of electrons to an electromagnetic 
field. 



The question now is (Pereira and Castro Neto 20091: 
How are the transport properties of Dirac electrons 
changed when transversing a region of strained graphene? 
As usual, the answer to this question is obtained by com- 
puting the transmission of the device by matching the 
wave functions from the left and right leads to those of 
the central region, at the positions x — Q and x — L. 
From the matching conditions we compute the total scat- 



tering matrix of the system (Nazarov and Blanter 2009 1, 



relating the incoming and outgoing waves. The scatter- 
ing matrix S is obtained easily from the total transfer 



18 



matrix of the structure Tg using the same formalism in- 



troduced in Sec. III.C In this case, the transfer matrix 
is given by 



T. 



1 
D 



u V 

V u 



piLqai 



ue 



'iLqj: 
iLq^ 



(29) 



-ie 



with ky, kx 



and qx 



where D = 4 cos cos 0, u = e^ + 
ianO = ky/kx, tan0 — [ky — 5)/qx, 
the transverse momentum, the longitudinal momentum 
in the leads, and the longitudinal momentum in the de- 
vice, respectively, and S = At/(vph). Finally, the energy 

in the leads has the form e = ^/fc^ -t- k^, and in the cen- 
tral region e — \J q^ + (ky — (5)^. The last equation shows 
that the effect of strain is to shift the position of the 
Dirac point in the Brillouin zone, a crucial effect on the 
explanation of the following results. Equation (29 1 was 



derived for energies in the continuum (no bound states). 
However, a fundamental property of the scattering ma- 
trix (or the Tg matrix for this purpose) is that bound 
states can also be obtained from the form derived for the 
scattering states, by looking at the poles of the S ma- 
trix. Since the 5— matrix is obtained from the inverse of 
the Tg matrix, its elements contain a factor which is the 
inverse of the determinant of T^. Additionally, the 5*1 1 
element of the 5— matrix (in this problem the 5— matrix 
is a 2 X 2 matrix, since we are working on the propa- 
gating mode base) gives the amplitude of transmission 
across the strained region, its value being 



5ii 



4 cos d cos d 



cos(L(7a;)(|u|2 - |i;|2) - isin(Lqj;)(|up - jup) 



(30) 

From Eq. pO| ) , we see that there are energies of perfect 
transmission, when Lqx = riTr, with n = 1, 2, . . .. On the 
other hand, considering the case fc^ = 0, corresponding 
to normal incidence on the boundary, we obtain for the 
transmission 



r = |5iip = 



(52cos2(Lga;) ' 



(31) 



which is a number smaller than one, meaning there is 
no Klein tunneling through strained graphene. At reso- 
nances, Lqx — nn, the transmission (31) is one, and for 



energies in the range 5/2 < e < 5, it decreases exponen- 
tially with L, since qx becomes imaginary. In general, 
for angles (between k and y— axis) satisfying the condi- 
tion Of > arccos(— 1 -I- 6/e) the transmission, as com- 
puted from Eq. ( pO| , is shown to be strongly suppressed 
for large L, an effect termed: transport gap formation 



(Pereira and Castro Neto 20091. True energy gaps can 



also be created in graphene by choosing an appropriate 



geometry of strain (Guinea et al. 20101. The reason be 



hind this transport gap mechanism is easy to understand 
given the discussion in Sec. ITTj in the strained region the 



electrons' wave function are no longer eigenstates of the 
helicity operator h. Therefore, since the helicity is not 
a constant of motion in this problem, Klein tunneling is 

lost and backscattering is allowed. 
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Figure 11 (Color online) Resonances, bound states, and sur- 
face states. Top panel: The parameters are L = 100 and 
6 = 0.1. In the resonances region, e > ky, the energy curves 
corresponds to the momenta q^ = mr/L. The horizontal axis 
refer to the transverse momentum ky, and the vertical one 
to the dimensionless energy e. Bottom panel: Representa- 
tion in the plane energy versus ky of the five types of states 
appearing in this problem. The two Dirac cones, that of the 
contacts and that of the strained region, are represented for 
kx = and q^ — 0, respectively. The off set of the apex of the 
two cones is S. All states with k^ j^ and Qx ^ lie above 
the two respective cones. 

As mentioned, the denominator of the Su element of 
the iS— matrix is all that it is needed to look for bound 
states in this system. We can imagine two different 
types of bound states: those decaying exponentially in 
the leads, but propagating inside the strained region and 
those decaying exponentially in the three regions. The 
latter states are edges states living at the boundaries be- 
tween the leads and the strained region. Looking only 
at the poles of Su we can find that both types of bound 



states exist (Pereira and Castro Neto 20091. The rich 



ness of states in strained graphene is due to the breaking 
of the chiral symmetry. In Fig. [TT] we have represented 
the energies at which the transmission is unity (called res- 
onances), the energies of the bound states, and the energy 
of the edge or surface states. The different types of states 
are a result of the different shifts of the two Dirac cones 



over the Brillouin zone (Hasegawa et al. 2006 Montam 
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baux et al. 2009 1 , induced by the strain and governing 
the vector potential of Eq. 28 the two apexes are shift 



by S due to strain, as seen above. 



It was shown (Pereira and Castro Neto 2009)that in 



the problem under study there are in total five types of 
states: (1) scattering states, (2) band states (states local- 
ized in the junction along the a;— direction) propagating 
along the y direction (3) localized states at the bound- 
ary of the junction, (4) filtered states, that is, scattering 
states decaying exponentially inside the junction for cer- 
tain values of the incoming angle Of, and (5) and states 
such that the transmission occurs via evanescent waves 
for any orientation of the incoming momentum - it is said 
that all states are filtered. The regions in the energy ver- 
sus ky plane where these type of states appear are shown 
in the bottom panel of Fig. [TT] 

In conclusion, the example discussed shows that there 
is plenty of room at the bottom of strained graphene for a 
whole new sub-field of graphene research: that of strain- 
based transport engineering or straintronics. 



V. QUANTUM CORRECTIONS TO THE DRUDE 
CONDUCTIVITY 

Before discussing the quantum corrections to the 
Drude conductivity, we introduce some key concepts on 
electronic transport. Consider first the elementary trans- 
port theory in a normal metal, where electrons have an 
effective mass m* . The velocity of the electrons at the 
Fermi surface is given by vp = hkp/m* , where kp is the 
Fermi wave number whose value depends on the density 
of electrons. Impurities in an otherwise perfect crystal 
occasionally deflect electrons from free propagation, lead- 
ing to the appearance of a mean free path I ~ the mean 
distance to the next collision. Since the dominant contri- 
bution to transport comes from electrons having veloc- 
ity vp, we can introduce a phenomenological parameter, 
the relaxation time t, defined as rvp = £; the elemen- 
tary theory of transport then shows that the electronic 
conductivity of the metal reads ( Ziman 19791 (no spin 



or valley degeneracies included): (Tq = frrvr/m* . This 
is a purely classical result, known as Drude 's formula, 
and which assumes that, after each collision, the electron 
looses memory of its previous linear momentum state. 
The calculation of r is usually obtained from Fermi's 
golden rule. The above description makes sense when 
(- ':^ \p ~ 2tt /kp (Ziman 19791. If we now repeat the 



same analysis for graphene, we obtain (Adam et al. 2008 



20071 |Basko| [2008| [Peres et al.\ |2007bt [Stauber et al 
2007): (To = 2e'^TVpkp/h, where we have now included 



the contributions of both spin and valley degeneracies. 

It is also possible to view the conductivity problem as 
a random walk. In this case, the conductivity is related 
to the diffusion constant D thro ught Einstein's relation 
(Chakravarty and Schimd 1986j ): ao — e^Dp{ep) (here 



again no spin or valley degeneracies included), with p{ep) 
the density of states per unit area, and the units of D 
are those of area per time, in any spatial dimension. The 
diffusion constant, for order of magnitude estimates, can 
be taken as Z? ~ vp£. 



A. Weak localization in a normal metal 

Weak localization is a correction to the classical con- 
ductivity of a disordered metal due to quantum interfer- 
ence, and originates in the quantum mechanical super- 
position principle. Electrons propagating in metals are 
subjected to a number of scattering mechanisms that give 
rise to a number of characteristic times. 

The relaxation time r due to elastic collisions with 
static impurities is assumed to be the smallest scattering 
time and describes a reversible process. Other scattering 
mechanisms are irreversible in nature and lead to either 
the loss of phase coherence or energy relaxation; for in- 
stance, those caused by electron-electron and electron- 
phonon interactions (excluding interactions with mag- 
netic impurities). Contrary to r, the phase relaxation 
(or dephasing) time T<p is temperature dependent and at 
low temperatures is mainly due to electron-electron inter- 
actions. In the presence of a magnetic field, a new time 
scale tb appears, which is of the order of tb ^ f^%/D, 
by simple dimensional analysis arguments, and where 
^s — fi/cB is the magnetic length. We now proceed 
to the discussion of the quantum interference effects us- 



ing an intuitive approach (Abrikosov 1988 Beenakker 



and Houten 19911, rather than a formal one, as seems 



appropriate in the context of this Colloquium. 



A type T 




Figure 12 (Color online) Pictorial representation for two 
types of scattering processes, whose physical interpretation 
is given in the text. The presence a finite magnetic field is 
represented by the fiux $ piercing the area defined by the two 
time reversed trajectories. 

Imagine an electron traveling from position AtoB, as 
shown in Fig. 12 and we denote by Oie"^' the probability 
amplitude for the electron to travel from A to B, along 
trajectory i. Since there are many indistinguishable tra- 
jectories, the total probability of traveling from A to B 
is 



Pi{A^B) 



E 



i(pi 



(32) 



We now show that quantum interference effects are much 
more important for what we call trajectories of type II, 
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in Fig. [T2j in which the initial and final points coincide 
iA = B). 
For type I trajectories {A ^ B), we have 



Pi{A^B) = 



Z^ 



11° 



jLljC 



.-0J-) 



(33) 



Since the phases for different trajectories of type I are 
uncorrelated, we assume that the second term averages 
to zero, leaving us with the classical result, in which the 
probability to go from A to i? is just the sum of the 
probabilities over all possible trajectories, that is. 



Pi{A -^B)=J2 W^f = Pi (A ^ B) 



(34) 



However, for trajectories of type II (the same initial and 
final points), in the presence of time reversal symmetry, 
the situation is quite different. In fact, time reversed 
trajectories (going round the loop in clockwise and anti- 
clockwise fashion), contribute to the sum of Eq. (33) 



with the same amplitude, in both modulus and phase. 
As a result, in the interference term. 






,) 



(35) 



when i and j denote time reversed trajectories, the phases 
cancel, even before any averaging; this term gives a con- 
tribution exactly equal to the first one, since for every tra- 
jectory there is a time reversed pair. This then amounts 
to a probability 



Pii{A ^ A) = 2^ |a,|' = 2P^j''\a ^ A) 



(36) 



This effect of quantum interference therefore enhances 
the probability of return, relative to the classical re- 
sult, decreasing diffusion and, therefore, the conductivity 



words we have 



(Abrikosov 1988 Beenakker and Houten 1991); in other 



(To <0, 



(37) 



where a^i stands for the conductivity of the metal con- 
sidering the enhanced backscattering effect, due to quan- 
tum interference. It is the reduction in a^^ji relatively to 
CTo that is known as weak localization. 

In the presence of a magnetic field B the relative phase 
of the electron's wave function, associated with the two 
time reversed trajectories of type II, has the value Scfi — 
47r$/(/>o as given by the Aharanov-Bohm effect, where 
3" is the magnetic flux piercing the area defined by the 
closed trajectory. Therefore, applying a magnetic field to 
the system suppresses the interference effect (because it 
changes the relative phase to a non-zero value) given by 



Eq. ( 36 1 , and the low-temperature conductivity of the 



metal increases when the field is turned on: or, in other 
words, we have 



Sa{B) _ awi{B) - a^ 



<7wl 



<7wl 



>0. 



(38) 



Using Eq. ( 38 1 we obtain experimental evidence of weak 



localization effects in the conductivity of a disordered 
metal. 



If the spin-orbit interaction (Hikami et al. 1980) can 



be ignored and there are no magnetic impurities in the 
metal, the times t, t^, and tb are the only relevant time 
scales, and they control the behavior of the low temper- 
ature conductivity. 

We now extend the previous analysis to graphene. 
Again, as expected, the chiral nature of the electrons (or, 
equivalently, their non-trivial Berry's phase) will play a 
major role. As before, we keep the discussion as elemen- 
tary as possible. 



B. Weak localization in graphene 

In graphene, the quasi-exact conservation of the chiral- 
ity and the existence of two valleys have profound effects 
in the low-temperature conductivity of the material. Be- 
low, we present the general picture of the quantum cor- 
rections in graphene, referring the interested reader to 
the literature for the subtleties appearing under a de- 



tailed analysis of this problem (Aleiner and Efetov 2006 



McCann etal] 2006| [Morpurgo and Gumeal 2006 Muc 






ciolo and Lewenkopf[|2010[ Suzuura and Ando[ 2002). 



As shown in Sec. V.A| in a normal metal the only 

elastic time is the relaxation time r. In graphene the 
situation is more complex. In order to understand the 
complexity of the situation, consider two different matrix 
elements of a potential created by a given impurity. We 
assume the potential to have the form 



V{r 



'Irl 



(39) 



The range of the potential depends on the value of tq: 
the larger ro, the larger the range. The effect of this 
potential on electrons within the same valley (denoted 
intra- valley scattering) is given by the matrix element 



(^+(fc')|V^(r-)|^+(fc)) = ^/(fc,fc')e 



-^rll^ 



(40) 



with /(fe,fc') = cos[6'(fe)/2 - 6l(fc')/2] and q = |fc - k'\. 
If we take ep = 0.5 eV one obtains q ^ Q.2/a, with 
a = -\/3ao. The function f{k,k') shows that the scatter- 
ing is not isotropic in momentum space, a consequence of 
the chiral nature of electrons in graphene. From /(fc, fc'), 
we also see that the scattering amplitude for backscatter- 
ing, /(fe, — fc), is zero, the fingerprint of Klein tunneling 
(recall Fig. [3]) for massless Dirac electrons (Beenakker 
20081. 
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If the potential also couples electronic momentum 
states from K and K' valleys (denoted inter-valley scat- 
tering) the matrix element, using wave functions from 
different valleys, reads 

(^+(fc')|^(r)|V+(fc)) = ^5(fc, k')e-Q''^o/i , (41) 

with gik,k') = isin[e'(fc)/2 - e{k')/2] and Q ~ \K - 
K'\ — 47r/(3a). In this case, backscattering is permitted, 
g{k, — fc) 7^ 0, since scattering couples states in the K 
and K' valleys, which have opposite chirality (recall Fig. 

i- 



Equation (41 1 shows that only for very short-range po- 
< a, does intervalley scattering have a signif- 
For long range potentials, only intra- 



tentials, r^ 

leant amplitude. 

valley scattering plays a role. 

From the above discussion it follows that, in the case 
of graphene, we need to define several elastic scattering 



Morpurgo and Guinea 2006): 



times (Aleiner and Efetov 2006 McCann et al. 2006 



1. Tiv, representing inter- valley scattering, whose scat- 
terers are very short-range potentials with range 



''o 



< a, such as some types of adatoms, adsorbed 
hydrocarbons, or vacancies. 

2. Tj, representing intra-valley scattering, whose scat- 
terers are long range potentials, such as ripples, 
dislocations, and charged scatterers. 

3. Tu,, representing also another contribution to intra- 
valley scattering. This scattering time has its ori- 
gin in the fact that chirality is not an exact sym- 
metry of Dirac fermions in graphene (due to trig- 
onal warping effects), therefore allowing for some 
amount of backscattering within the same valley. 
The importance of this scattering time grows as 
the Fermi energy increases. 

As long as the scattering potentials are long range, 
inter-valley scattering is negligible and back-scattering 
in graphene is absent, except from a small contribution 
from Tu,. The effect just described is, at the more fun- 
damental level, a consequence of the Berry's phase (of 
tt) acquired by massless Dirac electrons when they per- 
form a closed chiral orbit ([Berry and Mondragon |1987[ 



Mikitik and Sharlai 1999L The Berry's phase trans- 



forms the constructive interference, we described above 
for normal metals, into a destructive one, leading to weak 
anti-localization. Under these circumstances, the inter- 
ference effect seen in Sec. IV.AI for a normal metal cannot 
exist in graphene, forward scattering is enhanced and one 
should expect weak anti-localization effects to manifest 
themselves, that is, if r^ ^ t^, we have (ignoring t^u) 



(McCann et al. 2006 Suzuura and Ando 20021 



f^awl — (To > , 



(42) 



where aawi represents the enhancement (anti-localizing 
effect) of the conductivity over Drude's result ctq due to 
Klein tunneling. If short range scatterers are present, 
then inter-valley scattering plays a role, and since elec- 
trons in the K and K' valleys have opposite chiral- 
ity, back-scattering is present [see the function g(fc, fc')[, 
and we expect weak localization effects, according to Eq. 



37 detailed calculations confirmed this picture (McCann 



et al. 2006 Suzuura and Ando 2002). When the effect 



of a magnetic field is included, the general rigorous ex- 
pression for the weak localization corrections was derived 



by two groups independently (Aleiner and Efetov 2006 



McCann et al. 2006 1, and 5<t{B) can be either positive or 



negative depending on the relative values of the different 
scattering times, including r^. For a comprehensive dis- 
cussion of the interplay between the different scattering 
times and the quantum corrections to the conductivity, 
the interested reader is referred to the technical litera- 



McCann et al. 2006 1 



ture (Aleiner and Efetov 2006 Kechedzhi et al. 2007 
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Figure 13 (Color online) The quantity Sa{B) is defined as in 
Eq. (38 1. Left: Weak anti- localization behavior (the gate 
< 1 V, corresponding to an electron den- 



voltage used was Vg 

sity n < 7 X 10^" cm~^). Center: Weak localization behavior 
(the gate voltage used was Vg = 11 V, corresponding to an 
electron density n ~ 8 x 10^^ cm~^). Right: Dependence of 
Sa{B) on the electronic density (it grows from I to III), at 
T = 27 K. A fit of the data must use the rigorous formula 
derived in the literature (Aleiner and Efetov [2006 ' McCann 
" Wa l. , 2006). [Data from F. V. Tikhonenko et aL p"ikhoncnko 
et a/.||2008| |2009| ) ' 



In the left panel of Fig. [T3|we show, for small values of 
the magnetic field, a weak localization dip {5a{B) grows) 
followed by, above a certain field value B* , weak anti- 
localization behavior of the conductivity, since 6a{B) 
starts to decrease upon increasing the magnetic field over 
B* (no saturation of Sa{B) is measured upon increas- 
ing B, as in the weak localization case). In the central 
panel of Fig. |13| we show weak localization behavior 
in graphene, since the corrections 5<j{B) never decrease 
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upon increasing the magnetic field, and tend to satura- 
tion. Since the electronic density in the central panel of 
Fig. [13] is few times larger than that in the left one, it 
seems that the effect of short-range scatterers is more ef- 
fective at higher densities (Tikhonenko et al. 20091; at 



lower densities long-range scatterers dominate. Indeed, 
in the right panel of Fig. [13] we clearly see a crossover 
from weak anti-localization to weak localization as the 
electronic density increases from I to III, at a tempera- 
ture of 27 K; considering important screening effects of 
charged impurities at large electronic densities, such a re- 
sult sounds reasonable. It is a remarkable experimental 
fact that quantum interference effects just discussed can 
be observed in graphene at temperatures as high as ^-^200 



K (Tikhonenko et al. 20091 



As stated, we expect that in graphene the presence 
of different kinds of defects (in different concentrations) 
will control whether weak localization or weak anti- 
localization is observed. This depends on the relative 
value of the different elastic times introduced above and 
on the electron density. A detailed analysis of this point 
is essential for a correct interpretation of the data, and 



has been done in great detail ( McCann et al. 2006 Mor- 



purgo and Guinea 20061. The numerical values for the 



different scattering times can be obtained from the ex- 



perimental data (Tikhonenko et al. 2008). 



Finally, we note that the observation of quantum cor- 
rections to the conductivity in graphene seems to de- 
pend on the details of the fabrication process, which 
determines the amount of rippling introduced in the 
system ( jMorozov et al.\ [2006 1 . Routes for suppression 
of weak (anti-)localization effects have been considered 
(|Khveshchenko| 2006[) and this effect was experimentally 



observed as well (Morozov et al. 2006 Tikhonenko et al' 
[2008). As in the case of strain discussed in Sec. 
ripples are equivalent to effective gauge fields which break 
time reversal, leading to the suppression of weak localiza- 
tion effects, within each valley. In the system as a whole 
(both valleys considered) the full time reversal symmetry 
is preserved. 



VI. THE OPTICAL CONDUCTIVITY OF GRAPHENE IN 
THE INFRARED TO VISIBLE RANGE OF THE SPECTRUM 



In the ensuing sections we discuss the calculation of 
the percentage of light transmitted by a graphene mem- 
brane, when light shines from behind. This property is 
controlled by the optical conductivity a{uj) of the mate- 
rial. We analyze how and why the experimental behavior 
of a{uj) deviates from the predictions of the independent 
electron model. 



A. Graphene as a transparent membrane 

The calculation of light absorption by a given mate- 
rial is equivalent to the calculation of the optical con- 
ductivity. In general, such a calculation proceeds using 
Kubo's formula. In the case of graphene, it is possible 
to use Fermi's golden rule to obtain directly the fraction 
of absorbed light, which turns out to be a much simpler 
calculation than computing the optical conductivity first 



(Kravets et al. 20101. The central quantity to be com- 



puted is the transition rate of electrons excited from the 
valence band to the conduction one, as shown in Fig. [14] 
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Figure 14 (Color online) Pictorial description of the optical 
excitation of electrons in graphene. The absorption of a pho- 
ton can only induce vertical inter-band transitions. From left 
to right we have graphene doped with holes, neutral, and 
doped with electrons. 

In the presence of a vector potential A the Dirac 
Hamiltonian has the form 



Hk = vpcr ■ {p + eA) 



(43) 



We represent the electric field 'ds E = ~dA/dt and 
choose the polarization of the field along the a;— axis: 



A = xAo{e'' 



)/2. The term vpcr ■ eA will be 



taken as perturbation, and in the spirit of time depen- 
dent perturbation theory, only the exponential with neg- 
ative exponent is taken. The transitions induced by light 
absorption are now controlled by the ax matrix. Clearly 
the matrix element {ipx\o'x\'4'\) cannot contribute to the 
conductivity, since light cannot induce transitions within 
the same band, among states of equal momentum. The 
only non- vanishing contributing matrix element is there- 
fore {ipi\ax\ip-i) — —^vpeAosmO^k). The transition 
rate is then given by Fermi's golden rule: 



W^i,-i(fc) - 



2tt 



''Alsm^e{k)5{2vFkh~ujh). (44) 



The Dirac delta function in Eq. ( 44 1 enforces the condi 



tion that only electrons with energy uj/2 can be excited 
to the conduction band. The transitions we are referring 
to are shown in Fig. |14[ To obtain the contribution from 
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all states we have to integrate over the momentum and 
multiply the result by four (two for spin times two for 
valley). The calculations are elementary and the result 
for the total transition rate per unit area is 



1 _ e^Albj 



(45) 



If light of frequency w is shining upon a unit area of 
graphene, the amount of absorbed power per unit area is 
Wa = haj/r. The energy flux impinging on graphene is 
given by Wi — ceo-Eo/2, with Eq = Aqco. Therefore the 



fraction of transmitted light is ( Nair et al. 2008 1 



T = 1 



Wa 

W, 



= 1 



0.977, 



(46) 



with a ~ e^ /{AtteqIIc) the fine structure constant. The 
absorption of light is therefore independent of frequency 
and given only by universal constants. The high trans- 



mittance of graphene is shown Fig. 15 it is remarkable 
that a one atom thick membrane can be seen by the naked 
eye. 




Figure 15 (Color online) An optical image of an aperture 
partially covered with graphene and its bilayer (from left to 
right: air/graphene/bilayer), taken in a light transmission ex- 
periment (courtesy of A. K. Geim). 



It follows from the previous analysis that the transmis- 
sion at finite doping is given by 



Tiep) ~ {l-Tra)e{uj-2eF): 



where the Heaviside step function takes into account that 
absorption can only occur for frequencies larger than 
twice the Fermi energy, due to Pauli's principle. The 
result given by Eq. ( 46 ) is identical to the one given by 



a rigorous calculation based on Kubo's formula (Abergel 



eFaTI [20071 IGusynin e^ aT] [20071 [Kuzmenko et a?.||2008 



Peres et al.\ 2006 Peres and Stauber 20081. The reason 



why the perturbative calculation works so well is because 
the final answer is controlled by the small dimensionless 
parameter a. The reason why the transmission is con- 
trolled by the fine structure constant originates in the 



chiral nature of the electrons in graphene, a result exten- 



sible to few-layers graphene (Min and MacDonald 2009 1 



It is now a simple matter to include corrections to the 
linear spectrum of graphene in the formalism (the conical 
nature of the spectrum is valid for energies of the order 
of < 1 eV) . The addition of a next-nearest neighbor hop- 
ping term can also be included and treated within this 
formalism. The case of a next-nearest neighbor hopping 
is actually trivial, being proportional to the identity ma- 
trix its contribution is zero. In fact, it is proven in general 



( Stauber et al. 2008b I that the contribution of the next- 



nearest neighbor hopping only enters in the final result 
as a renormalization of the energy spectrum. 

As stated, the transmittance of light through graphene 
can be computed from a previous knowledge of the op- 
tical conductivity of the material. The transmittance is 
calculated from the solution of Fresnel's equations, read- 



ing ([ Abergel et aZ.[ [2007| [B lake et al] 2007 Pedersen 
[Stauber et aZ'l|2008b 



r=|l + a(c^)/(2c6o)r 



(48) 



where <7{uj) is the the optical conductivity (Falkovsky 



and Pershogubaj [2007[ falkovsky and Varlamov 



Gusynin and Sharapov[[2005||Gusyiiin et al.\ 2007[ 



2007 



2009 



Peres et al] [2006[ [Stauber et al] [2008b[ ) of graphene, 
given, at zero temperature and within the independent 
electron approximation, by (ep > 0) 



cr[uj) = aoO{ijjh — 2ep) 



4e, 



lao- 



.CTo \hui- 
TT \riuj - 



2e, 



2eF\ 



The quantity ag = 7re^/(2/i) is termed the ac universal 
conductivity of graphene. Inserting Eq. (49 I into Eq. 



(48 1 and taking cf = we obtain the result of Eq. (46 1 



Working the other way around, using Eq. (47 1 in Eq 



(48 1, we obtain for neutral graphene (t{uj) — aoi in accor- 



dance with Eq. ( 49 ) . The fact that for neutral graphene 



cr{uj) is given in terms of universal constants only, with 
no reference to any of the material parameters, is a rare 
result in condensed matter physics. If the intensity of 
light impinging on graphene is large, then non-linear cor- 
rections to T start to play a role, which is expected to 



lead to an increase of the transmittance (Mishchenko 



(47) 2009 Rosenstein et al. 2010). The non-hnear optical 



susceptibility coefficients of graphene have recently been 



measured using four-wave mixing (Hendry et al. 20101 



Due to its high transmittance, graphene can be used as 
transparent conductive electrodes in solar cells and liquid 



crystal devices (LC D) ( [Bae et aZ.[ [20101 [Blake et a;.[[2008 
Wang et aZ.[ [20081. 



B. The optical conductivity of neutral graphene 

We now discuss the experimental results for the opti- 
cal conductivity of neutral graphene and how those mea- 
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surements deviate from the independent electron model 
presented above. 



2006 1 calculation of the conductivity of graphene at zero 




0.8 1.2 1.6 2 

photon energy (eV) 



Figure 16 (Color online) Transmittance of graphene. For pho- 
ton energies below 1.2 eV, we plot data from reflectance mea- 
surements taken by K. Fai Mak ei al. (Mak et al. 20081; dif- 



ferent curves correspond to different graphene devices (cour- 
tesy of F. Kin Mak). For photon energies larger than 1.2 eV, 
we plot data from two different types of measurements (|Nair| 



et al. 



20081): narrow bandpass filters, with a full band width 
at half maximum of 10 nm, and standard spectroscopy mea- 
surements. The dashed line is the result of Eq. ( 46 1 . The solid 



(green) line, to the right of 1.2 eV, is a plot of T obtained from 
a calculation of the conductivity, including trigonal warping 



corrections ( Stauber et al. 2008b I 



In Fig. [16] we show measurements from two differ- 



ent groups ( |Mak et al] [2008} |Nair et al.\ |2008[ ) in two 
different regions of the spectrum. For photon energies 
E^ below 1.2 eV (in the near infrared region), the data 
were taken from reflectance measurements, at room tem- 
perature. The absorbance of the different samples is 
spectrally flat within a band of 10%. Roughly speak- 
ing, the transmittance follows the value given by Eq. 
( 46 1 . A more detailed analysis shows that the trans- 
mittance below E^ <0.5 eV increases over the universal 
value, whereas closer to E^ ^1-2 eV it decreases slightly 
from that value. According to Mak et al. (Mak et al. 



2008), and for energies E^ <0.5 eV, both temperature 



effects and some amount of variable extrinsic doping de- 
crease the conductivity and therefore produce an increase 
in the transmittance. In other words, if light absorp- 
tion decreases, then the transmittance increases. This 
effect is equivalent to a finite chemical potential (hav- 
ing effectively a Pauli's principle based blocking effect), 
which contributes to an increase in the transmittance. 
The temperature effect that they ( |Mak et al.\ |2008| used 
in their argument, was predicted to be of importance at 
room temperature and energies below 0.5 eV for undoped 



graphene (Peres and Stauber 2008), in agreement with 



the experimental measurements. Additionally, this set 
of measurements showed that the optical conductivity 
of graphene can also be affected by unavoidable doping 



temperature and zero doping showed that cr(w) departs 
from its universal value, a^ — ire^ /{2h), being strongly 
reduced at low frequencies, due to disorder. 

In conclusion, the deviation of the transmittance of 
graphene from the universal value predicted for Dirac 
fermions is a way of gaining insight on other electronic 



effects present in the material (Li et al. 2008b). 



Figure [16] also shows the transmittance of graphene in 
the photon energy range 1.2-3 eV. The measured value 
follows the prediction of Eq. (46), except at energies 
around 3 eV, where absorption increases. Since at ener- 
gies as large as 3 eV the electronic energy dispersion de- 
viates considerably from the Dirac cone approximation - 
an effect known as trigonal warping - a calculation tak- 
ing into account trigonal warping corrections to the band 



structure of graphene was performed (Nair et al. 2008 



Stauber et al. 



2008b), and the result is given by the solid 



line in Fig. 16 The calculation does predict an increase 
in light absorption at energies around 3 eV (an effect 
opposite to that discussed above for energies E-^ <0.5 
eV), but falls short on accounting for the magnitude of 
the effect. Within the independent electron model, the 
computed enhancement of the conductivity is essentially 
due to the increase in the density of states as the Van 
Hove singularity is approached, which is located at the 
energy of ^2.7 eV. Thus, transitions between states lo- 
cated at the Van Hove singularities of the valence and 
the conduction bands require photons of energy ~5.4 eV. 
Optical experiments using photons of about this energy 
have confirmed the strong enhancement of light absorp- 



tion due to the Van Hove singularities (Kravets et al. 



2010) 



Two possible additional causes for such an increase in 
the absorption come to mind: contamination of the sam- 
ple due to some organic residues (originated from the 
exfoliation process) or/and many-body effects (Herbut] 



et al.\ |2008[ [Mishchenko] [2008[ ). A recent calculation, 



however, showed that electron-electron interactions cor 
rect the transmission (reducing it) by only 0.03 — 0.04% 
( [Katsnelson"! [2008[ [Sheehy and SchmaHan[ [2007] [2009 ). 
On the other hand, in the regime of energies relevant 
for the visible range of the spectrum, the fine structure 
constant of graphene ag is a number of order 1, a fact 
casting reasonable doubts on the validity of any pertur- 
bative calculation. 



A recent ah initio calculation (Yang et al. 2009) used 



and intra-band scattering. Indeed, a SCBA (Peres et al. 



the Kohn-Sham eigenvalues and eigenvectors to compute 
the optical conductivity of graphene from infrared to vis- 
ible frequencies. This type of mean-field calculations 
includes electron-electron interactions by means of the 
exchange and correlation approximations. The optical 
conductivity obtained does not fit exactly the data of 
Fig. |16[ since it predicts an absorption higher than what 
is measured experimentally, in the full frequency range. 
Nevertheless, the proposed excitonic effects (included via 
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the solution of the Bethe-Salpeter equation) were shown 
to account well for the experimental deviations (to lower 
energies) in the enhancement in light absorption for pho- 
ton energies associated with transitions between states lo- 
cated at the Van Hove singularities of the valence and the 



et al. 2010 



In Fig. 17 we show one set of the mea- 



conduction bands (Kravets et al. 2010) relatively to the 



predictions made by the non-interacting electrons the- 
ory. Indeed, the elementary theory predicts an intense 
absorption peak at about 5.2 eV, whereas the measured 
data shows a red shift of the peak to 4.6 eV. The discrep- 
ancy can be solved by considering the mutual attraction 
of the electron-hole pair, created at the two Van Hove 
singularities, when a photon of the right frequency is ab- 
sorbed. 

The energy red-shift mentioned is easy to understand 
from the point of view of the mutual attraction of the 
electron-hole pair (forming an exciton) : since the electron 
and hole have opposite charges, the effective energy seen 
by the photon being absorbed is the non-interacting value 
minus a positive correction coming from the electrostatic 
attraction between the electron and the hole created by 
the absorption of the photon. 



C. The optical conductivity of gated graphene 



sured data (solid blue curve), together with calculations 
using two models: (i) a model where the independent 
electron theory is supplemented with the effect of dis- 
order (computed with the two models discussed in Sec. 
IV. A I and (ii) a calculation including excitonic effects on 
top of model (i). It is clear that the first one, considering 
only disorder, does not account for the five deviations 
observed in the data, relatively to the independent elec- 
tron model; it partially accounts for the enhancement 
of the conductivity below 2/x. The additional optical 
response observed in the data, in this frequency range. 



must come from intra-band scattering processes ( Grushin 



et al. 2009j), not included in model (i). Points (ii)-(w) 
are all accounted for by including excitonic effects, that 
is, model (ii). There is, however, the need to use in the 
model a higher temperature, T — 120 K, than that mea- 
sured experimentally by the finger of the cryostat, T — 54 
K. This can be justified by recalling that graphene is a 
system presenting small inhomogeneities of the electronic 
density. This fact effectively smears the chemical poten- 
tial and such an effect can be accounted for by consider- 
ing an effective higher temperature, as was shown to be 
the case in the interpretation of the optical response of 



graphene's bilayer (Kuzmenko et al. 20091. 



Measurements of the optical conductivity of gated 
graphene, in the far infrared region of the spectrum ( Li 



et al. 2008b I, also show strong deviations from the sim- 



ple theory given above, and expressed in concise form by 
Eq. ^ 
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The deviations seen in the data (Li et al. 2008b I, for 



all values of the gate-voltage considered in the experi- 
ment (ranging from Vg =10 to 71 V), are of five different 
types: (i) finite absorption below 2/x, which is due to both 
inter-band and intra-band elastic and inelastic scattering 
processes; (ii) broadening of the absorption edge around 
the energy threshold 2/x; (in) an enhancement of the con- 
ductivity above the universal value ctq hi the energy range 
between 2/i and 2fj, + E* , where E* is a characteristic en- 
ergy scale, (iv) a reduction of the conductivity bellow ctq, 
at energies above E* , with the conductivity as a func- 
tion of frequency having a positive curvature, and (v) 
the imaginary part of the conductivity is larger than the 
value predicted by the non-interacting model for energies 
hv > 2/i. 

Additionally, we point out that the optical conductiv- 



ity curves measured experimentally ( Li et al. 2008b ) col- 



lapse on top of each other when re-plotted as function 
of Lu/ fi, implying that the mechanism causing deviations 
from the non-interacting approximation must be intrin- 
sic. 

To explain all measured deviations, it is necessary to 
take into account disorder, temperature, electronic den- 
sity inhomogeneities (Zhang et al. 20091, and electron- 
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Figure 17 (Color online) Optical conductivity of gated 
graphene on top of silicon oxide, in units of ao. The data 
is given by the solid curve. The dashed curve is a calcula- 
tion taking disorder into account, and the dot-dashed curve 
is a calculation including disorder and electron-electron in- 
teractions. The experimental curves refer to a gate voltage 
of 28 V, which corresponds a Fermi energy of ^ ~ 0.18 eV. 
The calculation took an effective temperature of 120 K. In 
the left (right) panel we have the real(imaginary) part of the 
conductivity. In the inset, the measured dc conductivity of 
the device (solid curve) depicted together with a fit of a{eF) 
using Eqs. ( 20 1 and ( 25 I. The fit of the dc conductivity fixes 



electron interaction effects, of excitonic nature (Peres 



the concentration of impurities to the values of rii = 2 x 10 
cm~^, for resonant scatterers, and nt — 1 x 10^^ cm~^, for 
charged ones (dashed curve). (Data from Z. Li et al. ( [Li| 
\et al.\^008h\, courtesy of Zhiqiang Li.) 
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In the same set of measurements (Li et al. 2008b I, an 
increase in the Fermi velocity over the value 1.1 x 10^ 
m/s was measured, upon diminishing the gate voltage. 
The renormalized Fermi velocity, due to exchange, was 
predicted to be (including the Thomas-Fermi screening 
of the Coulomb potential) 



VF = "ybarc + 



47reoed^ 



In 



2qc 



(50) 



where Vbare is the bare Fermi velocity, and therefore not 
an observable parameter, Qc ^ I A~^ a cutoff momentum, 
and we have used the fact that the Thomas-Fermi screen- 
ing momentum is for graphene on silicon oxide given by 
(7tf — 2fcF [we note that Eq. ( 50 1 is an extension, tak- 



ing screening into account, of previous results (Gonzalez 
efal] 1 1996^, Hwang et a;.|[2007l|Polini ef'ZjpOTTl )]. The 



result ( 50 1 does show that vp increases upon decreasing 
kp, but the formula fails to fit the data (Li et al. 2008b) 



over the measured gate voltage range, especially for larger 
Vg. Alternatively, phonons also seem to partially account 
for the Fermi velocity renormalization seen in the exper- 



iments (Gusynin et al.\ 



et al. 2008a I 



2009^ Peres et al. 



2008 Stauber 



The renormalization of the Fermi velocity suggests 
that contributions from many-body effects can be impor- 
tant, but an experiment with a higher degree of accuracy 
should be repeated, since cyclotron mass measurements 
( [Novoselov et al. 2005a I do not see deviations of the 
Fermi velocity from vp — 1.1 x 10^ m/s upon varying the 
electronic density. 

Finally, we note that, as in the case of neutral graphene 
discussed in Sec. |VI.B| the electron-hole pair mutual- 
attraction also induces a shift in the Fermi edge at w = 2/i 
toward lower energies (red shift effect), since, as before, 
it renormalizes the non-interacting energies by adding a 
negative number. 

In conclusion, the optical response of graphene is a 
property of the material where different kinds of quantum 
effects seem to play an important role, all on equal foot- 
ing. Exploring graphene for nanophotonic devices surely 
requires a detailed understanding of its optical proper- 
ties. 



VII. CONCLUSIONS 

We have given an introductory review on the transport 
properties of graphene, touching the relevant aspects of 
this topic. A survey of the literature was given along with 
a discussion of elementary models, which, despite their 
simplicity, are in good quantitative agreement with many 
features of the data. Many of these models were discussed 
in greater detail in the literature (to which due credit is 
given), but using more formal methods. We hope that 
the topics discussed in this Colloquium may contribute to 



a wider dissemination of the physics of graphene among 
the non-specialized audience. 
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